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THE  APPLICATION  OF  A  COMPUTER  METHOD  TO  TIE  ANALYSIS  OF 
A  GRIDDED  CIRCULAR  CURVED  FRAME 


By  Shyi-Jiun  Leou 

SYNOPSIS 

The  analysis  of  the  gridded  circular  curved  frame  can  be  simplified 
by  using  the  computer  method.  Formulas  are  supplied  by  applying 
Castigliano's  theorem.  They  show  the  relationship  between  the  end  moments, 
end  torques  and  the  corresponding  rotations  and  displacements  of  the 
member  end.  A  proposed  matrix  method  associated  with  computer  analysis  is 
used  in  the  procedure  to  solve  the  problem.  The  use  of  this  method  is 
illustrated  by  a  detailed  numerical  example. 


INTRODUCTION 

To  meet  the  requirements  of  the  modern  trend  in  architecture  and  the 
clearance  and  space  of  some  special  structures,  it  seems  that  there  will 
be  an  increasing  use  of  gridded  curved  frameworks. 

Because  of  the  relative  influence  of  deformation  between  curved 
girders  and  cross  beams  throughout  the  whole  framework,  the  analysis  of 
the  gridded  curved  frame  is  a  difficult  task. 

In  this  report,  the  computer  method  will  be  used  to  analyze  such  a 
framework.   It  allows  the  greatest  possible  simplification  in  the  computa- 
tion required.  The  member  can  have  any  curvature  in  the  plan,  including 
segmental,  but  particular  attention  will  be  paid  to  the  common  circular- 
arc  shape.  Formulas  are  presented  for  frames  of  constant  cross  section 
that  give  the  relation  between  the  end  moments,  torques,  and  shears,  and 
the  corresponding  rotations  and  displacements  of  the  end  sections.  These 
relations  have  been  derived  from  the  expression  for  the  total  strain  energy 
due  to  bending  and  torsion  by  applying  Castigliano's  theorem.  From  these 
relationships  the  stiffness  matrix  of  individual  members  can  be  written  and 
the  stiffness  matrix  for  the  whole  structure  follows.  Through  computation 
with  the  electronic  computer  the  end  forces  can  be  obtained.  These  end 
forces  are  combined  with  the  fixed-end  forces  to  obtain  final  end  forces. 


STIFFNESS  MATRIX  OF  AN  ARC-SEGMENT 

A  beam  curved  in  plan  with  end  "a"  free,  and  end  "b"  fixed,  is  shown 
in  Fig.  1.  Assume  that  the  loading  and  translation  are  all  perpendicular 
to  the  plane  of  the  arc-segment. 


Fig.  1.  A  beam  curved  in  plan  with  one  fixed  and  one  free  end 

At  end  "a"  (Fig.  1),  a  moment  M^,  a  torque  T  ,  and  a  vertical  force 

F 
F  ,  are  acting  on  the  beam  as  indicated  by  the  vectors. 

The  properties  of  segment  ab  can  be  described  by  the  partitioned 

matrix  equation, 
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D  indicates  displacement  in  the  direction  perpendicular  to  the  plane 
at  end  "a".  Rra  and  R,  indicate  rotation  about  radius  and  tangential  axes 
at  end  "a"  respectively. 

Then  the  K's  in  the  stiffness  matrix  can  be  obtained  as  follows: 
At  any  section  an  angular  distance  9  from  "a"  the  forces  can  be  expressed 
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M„  =  FFrsin9  +  Hsa 
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TQ  •  Fr(l-cos9)   +  lysine  +  TFcos9. 
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A  sign  convention  that  can  be  applied  to  the  end  of  any  span  is  shown 
in  Fig.  2. 

The  moments  are  positive  when,  looking  outwards  from  the  center  of 
curvature,  they  produce  a  clockwise  rotation  of  the  section. 

The  torques  are  positive  when,  looking  along  the  tangent  to  the  beam 
in  a  counterclockwise  sense,  they  cause  a  clockwise  rotation  of  the 
section. 

The  vertical  forces  are  positive  when  directed  upwards. 


Fig.  2.  A  sign  convention  that  can  be  applied  to  the 
end  of  any  span 

The  total  strain  energy  stored  in  the  beam  due  to  the  acting  forces, 

neglecting  the  effect  of  vertical  shear,  will  be 
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EI  and  GJ  are  the  flexible  rigidity  and  torsional  rigidity 
respectively. 

The  vertioal  displacement  D   at  end  "a"  in  the  direction  of  T    is, 


by  Castigliano's  theorem,  I) 
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Then 
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Integrate  uid  get 
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p 
+  -igj-  r  [a  3in24>    +  2K(2  -   2oos^>  -  sin2o>  )] 

2 
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The  rotation  R       at  "a"   in  the  direction  of  vF  is  R       ■  -—-=-   . 
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The  rotation  R.   at  end  "a"  in  the  direction  of  /  is  1   -  -—-*  . 
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If  the  following  substitutions  are  made  in  the  above  equations, 
a11   =  4E1  [2$~  Bin24>  +  K(64>+  sin2£-  8ain<j>)]    . 
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r  21 

2 
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Equation  (4),  (5)  and  (6)  can  now  be  expressed  in  the  form  of  a  matrix: 
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Equation  (9)  shows  that  (A)"  is  a  transformation  matrix  between 

displacements  and  forces  at  end  "a",  but  this  is  what  the  matrix  K*.  in 

ab 

Eq.  (l)  stands  for.  So  we  get 


*ab  -  (A)"1- 


If  we  rewrite  Eq.  (2)  in  the  matrix  form,  we  get 
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Then  the  stresses  at  end  "b"  can  be  obtained  by  substituting  the  values 
of  F  ,  W    and  T  into  Eq.  (10)  and  by  changing  the  angle  6  intod>,  or 
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Thus 
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To  determine  the  values  of  K,   and  K  ,  we  keep  end  "a"  of  the  curved 
beam  in  Pig.  1.  fixed,  and  end  "b"  free.  Then  a  set  of  stress  equations 
for  the  section  at  an  angular  distance  9  from  end  "b"  can  be  obtained, 
and  are 
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Similarly,  following  the  steps  from  Eq.  (3)  to  Eq.  (7),  we  get 
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Each  element  of  the  B  matrix  in  Eq.  (14)  is  equal  to  each  corresponding 

element  of  A  matrix  in  Eq.  (7),  but  the  signs  are  not  all  the  same.  The 

elements  b .,.,  b  ,  b„„,  b   and  b,,  have  the  same  sign  as  the  element  a.., 

a1V  a22'  a51  anci  a5V  l)ut  b12'  b21'  *2^  anli  N2  haVe  the  °PPos:i-'';e  ai-Sn   to 
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Comparing  Eq.  (15)  with  Eq.  (1),  we  know  that 

*l  "  (B)"1 " <16> 

The  stresses  at  end  "a"  can  be  obtained  by  substituting  FT ,  (C  ,  T? 

into  Eq.  (13)  and  by  changing  the  angle  9  into£. 


So  we 
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As  an  example  the  stiffness  matrix  of  a  whole  structure  is  presented 
follows: 


Fig.  3.  Plan  view  of  a  frame  with  curved  and  straight  members 
Figure  3  indicates  a  frame  with  curved  and  straight  members  meeting 
at  joints  1,  2,  ...  ,  8.  Numbers  in  parentheses  in  Fig.  3  indicate  the 
numbering  system  for  the  bars.  The  relationship  between  member  end  forces 
and  elastic  member  end  displacements  can  be  indicated  by  Eq.(l9). 

{pj  '   [K]  M (19) 

Vector  P  in  Eq.  (19)  represents  member  end  forces.  Its  transpose  is 

JP)T  .  ,p1,p1  p4  p2  p2  P5  P?  P3  p6  p6  4  p7  p7  p5  p8  p8  ,     » 

V]  [VVVVV  3'  3'  4'  5'  e'V  6**7'  7'  7'  8j  *  (20' 

Symbol  P^    indicates  the  forces  at  end  l(the  subscript)  of  member  (l) 
(the  superscript). 

For  each  value  of  P  on  the  right  side  of  Eq.  (20)  there  will  be  three 
components.  For  example, 


-(21) 


On  the  right  side  of  Eq.(2l)  components  indicate  vertical  force, 
moment  and  torque  at  end  3  of  member  (2). 


11 


Vector  D  in  Eq.(l9)  represents  member  end  displacements.   Its  transpose 


la 


W  '  "   £>^^444^4Dt^D64DrDrD|j (22) 

For  each  value  of  D  on  the  right  side  of  Eq.(22)  there  will  be  three 
components.  For  example, 
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■(23) 


On  the  right  side  of  Eq.(23),  components  indicate  vertical  displace- 
ment, rotation  about  the  radius  and  rotation  about  the  tangent,  respectively, 
at  end  2  of  member  (4). 

The  stiffness  matrix  K  in  Eq.(l9)  is  an  end  displacement-end  force 
transformation  matrix.  It  is  shown  in  Eq.(24)  in  which  for  example  Kt 
is  the  stiffness  matrix  for  member  (2) (the  superscript),  and  it  relates 
the  displacements  at  end  2(the  first  subscript)  and  the  forces,  which  are 
induced  by  the  displacements  at  end  2,  at  end  }(the  second  lower  subscript). 


12 


r--      go 

CD  03  CO  CO 


IT  "IT 


"it 


"it 


*it 


MO  \D  MD  MD 


CM 


MD 


*3-  cm                                ^J-  CM  -rt-  CNJ 

Q  •       W 1*3 

«-  CM 

i-   C\J       i-   CM   i—  CM 

«         w       ^ 

*-  CM 

^q^-    *V"  *V «... 


r  CM  CM   CM   CM   N*,,ir\N'\KVN-\Kv*'<a-MD   IPvMDMO   ^J-MD   f— MD   f^-  f-  IT\  f—  CD  r- 

c^fx.c^pHfLtcmcMo^eMfi.cup^pL, 


13 


FIXED-END  STRESSES  DUE  TO  UNIFORM  LOAU 

The  end  stresses  for  a  beam  curved  in  plan  with  both  ends  fixed  and 
loaded  uniformly  with  q  per  unit  length  of  beam  will  now  be  determined. 


Fig.  4.  Beam  curved  in  plan  with  both  ends  fixed. 


From  Fig.  4,  the  vertical  reaction  F   at  end  "a"  is  equal  to 
a      2   ' 


and  the  distance  OG  from  the  center  of  curvature  to  the  center  of  gravity 
of  arc  ac  is 
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At  a  distance  r9  from  end   "a", 
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2                  "5" 
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Applying  a  unit  moment  M'  ■  1  at  "a".  Then,  on  the  section  at  an 
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angular  distance  9  from  end  "a",  Mi  -  cos9  and  T„  -  sin9. 
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As  the  team  is  fixed  at  end  "a", 

or 
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+     J    -gj-   [^ine  +  T^cos9  +  -3£.£(i_0ose)  -  qr2(9-sin9)J    sin9  de  =  0. 

Integrate  and  get 
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If  a  unit  torque  T£»  =  1  is  applied  at  end  "a",  on  the  section  at  an 
angular  distance  9  from  end  "a",  M£'-  sine  and  T£'-  cos9.  As  the  beam  is 

fixed  at  end  "a", 


r  J° 

or 
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Integrate  and  get 

M^(K-l)sin2<£  +  /[$(K+1)   +   (K-1  )sin*cos£]   -  qr2   [-l-£2(K+l) 

+  -5-^(K-l)sin^cos^     -  (K-l)ain2d5     +  K£sind>    -   2(K+l)(l-cos£)l   -  0. 


(28) 

From  Eq.  (27)  and  (28),  the  values  of  M^  and  T  can  be  obtained  for  any 

a     a 

known  values  of  the  uniform  loading  q  and  the  angle  4).   It  will  be 
convenient  to  make  the  following  substitutions  in  the  above  equations: 
'A  =  $(K+1), 
B  =  (K-I)sin2£, 
C  m   (K-l)sin$cos£, 

D  -4>(K+1)   +  K<j<cos<f>  +   1)   -  -1-4>(K-I)sin2^>,     (29) 

-  2(K+l)sin<£-  (K-1 )sin£cos£, 
_E  =  -j-4>2(k+1)  +  -j-^(K-l)sin^cos^  -  (K-l)sin24> 
+  K^sin4>  -  2(K+l)(l-cos£). 
Eq,  (27)  and  (28)  can  now  be  expressed  in  the  forms 

(A-C)^  +  BT1  +  qr2D  =  0,  and  (30) 

nf  +   (A+C)'/  -  qr2E  =  0.  (31) 

F      F 
Solving  for  M  and  T  gives 


„F  qr    [D(A+C)   +  BE] 

a  -  -  -?-:-c-2-:-bS— 

_F  3r2[E[A-C)   +  BD] 


and 


-(32) 

-(33) 
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FIXED-EN1)  STRESSES  DUE  TO  A  CONCENTRATED  LOAD 

The  end  stresses  for  a  beam  curved  in  plan  when  both  ends  are  fixed 
and  loaded  with  a  concentrated  load  located  at  an  angular  distance  oC 
from  end  "a"  will  now  be  determined. 

In  Fig.  5,  let  the  displacements  at  end  "a"  due  to  end  stresses  only 


ID 


be  (lW,l  -  <  D   L  ,  and  let  the  displacements  at  end  "a"  due  to  the  ac 


ting 


taJ 


concentrated  load  P  only  be  JDp]  . 

Then  according  to  the  fixed-end  condition  at  "a",  we  get 


{DFNS}   +  {Dp}   -°" 
By  Eq.(8),  we  get 


(A) 


Fi 


'•  -hf 


-(34) 


-(35) 


a 


(A)-1{Dp| 


-(36) 
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Fig.  5.  Plan  view  of  a  curved  beam  with  both  ends  fixed  with 
a  concentrated  load  P  acting  downward.  Oidioi 


From  Fig.  5>  we  find  that  the  total  stresses  on  the  section  at  an 
angular  distance  9  from  end  "a",  when  O^O^oc  ,  will  be 


9 


_F 


VL  •  Mcos9  -  Tsin9  +  Frsin9, 

w    a       a  a     ' 

TQ  ■  Msin9  +  T  cos9  +  Fr(l-cos9), 

9    a       a  a  v      ' 


-(37) 


From  Fig.  6,  the  total  stresses  on  the  section  at  an  angular  distance 
9  from  end  "a",  when  tx±9±$   ,  will  be 


V   =  F        -  P 
9    a 


M.   -  Mcos9  -  Tsin9  +  F'rsin9  -  Prsin(9  -o<), 
9         a  a  a  \  /i 

TQ  -  JCsin9  +  Tcos9  +  Fr(l-cos9)   -  Prfl   -  cos(8-<x)] 


-(58) 
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Fig.   6.     Plan  view  of  a  curved  beam  with  both  ends  fixed  with 
a  concentrated  load  P  acting  downward.    «iQi  j> 


If  the  effect  of  shear  is  neglected,   the  total  strain  energy  is 


U  -]0-2E?-  rd9  +/o-5§r  rd9  +Ja-5|r  rd9  +£  -^  rd9   . (59) 


Applying  a  unit  force  F'   .1   at   "a".     Then,   on  the  section  at  an 
angular  distance  6  from  the  end   "a",   Ml   ■  rsin©  and  T'    -  r(l-cos9)  when 
0£9i£.     As  the  beam  is  fixed  at  end   "a", 


-*H     =  D        ■   C"  M>    M     rd9      *  !X   T>    V      rd9 
~p       Dza  ")0  M9  VeI"  +J0    T9  VST" 


^v^-^v^--*. 
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J    "EI""   L^0039  "  Tasin9  +  Ursine]     sin9d9 

+     Jo~Sr     [^ainS  +  T^oose  +  P^r(l-cos9)]       (l-cos9)d9 

2 
+    i«  "EI"  K0039  -  Ta3^9  +  Far3in9  -  Prsin(9-«)]      sin9d9 


2 

r 

r 

E 
2 


j«  "ST     K8in9  +  Ta00s9  +  F^r(  1-oos9)   -  Pr  1-oos(8-«)J    (l-cos9)d9  -  0. 


Collecting  terms,   we  get 
2 


""Il""a[io^r3in  9  +  &  "  2Kroos9  +  Kroos29)d9 
r 

+  JK(rsin  9  +  Kr  -   2Kroos9  +  Kroos29)d9] 

2 

+  "Ei"  fr[  Jo(sin9cos9  +  Ksin9  -  Ksin9oo39)d9 
r 

+    JK(sin9oos9  +  Ksin9  -  Ksin9oos9)d9] 

2 
r         Fr    C,  2  2 

+  -gj"  Ta|_  J0("sin  9  +  Koos9  -  Koos  9)d9 
r 

+     JK(-sin  9  +  Kcos9  -  Koos29)d9] 
2 

■  "  _EI~[  J    (-Prsin9sin(9-«)   -  KPr  +  KPrcos  ( 9- «) 

r    " 

+  KPrcos9  -  KPrcos  ( 9- «)cos9d9]  .  (40) 

Apply  a  unit  moment  M^'  =1  at  "a".  Then,  on  the  section  at  an  angular 
distance  9  from  the  end  "a",  M£>  -cos9  and  T£>=sin9  when  0*9*^>. 
As  the  beam  is  fixed  at  end  "a", 

mf    "ra    J  o  9  Ver  +  J  0  T9  Vcr- 

a 
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f0*     r       r  F  F  Ft 

-b»—      M  cos9  -  T  sin9   +  F  rsinQ      cos9d9 
Jo   hi       L  a  a  a  J 

+  [    -*j--JMsin9  +  t'cos9  +  F  r(l-cos9)l    sin9d9 

+   f   -gj-- [M^cose  -  T^sin9  +  F*rsin9  -  Prsin(9-«)]    cos9d9 

+  J    -g-j-JM^3in9  +  T^cos9  +  F^r(l-cos9)  -  Pr  [l-cos(9-«)]  }sin9d9  .  0. 


Collecting  terras,   we  get 

ST      a 

r 


Fal  I     (rsin9cose  +  Krsin9  -  Kr3in9oos9)d9 

+    |     (rsin9oos9  +  Krsin9  -  Kr3in9cos9)d9] 
+  -gi~M^  [["(cos^  +  Ksin29)d9  +  [     (cos29  +  Ksin29)d9] 

f  r  c**  ($ 

a.   U  (-sine°°se  +  Ksin9oos9)d9  +  J   (-3in9oos9  +  Ksin9cos9)d9l 


"EI" 


"El"  J   [-Prsin(9-<x)cos9  -  KPrsin9  +  KPrcos(9-o()sin9]d9. 


-(41) 


Apply  a  unit  torque  T'"  =  1  at  end  "a".  Then,  on  the  section  at  an 
angular  distance  9  from  the  end  "a",  M'"  •  -sin9  and  T'"  -  cos9  when  0^9^. 
As  the  beam  is  fixed  at  end  "a", 

?qAF   "ta  )0"e     "BEI.    )„     B     ^9  GJ 
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j^-p.j—  [hToobO  -  T^sin9  +  F^rsine]   (-sin9)d9 
+  J    -gj-    [Masin9  +  T^ose  +  F^rO-cosS)]    cos9d9 

+  i«"EI""[  a  osS  "  Tasi"e  +  ^rsinS  -  Pr8in(6-«)]    (-sin9)d9 
r 

+  j«"  5J"{  I"Cain9  +  T^cose  +  Far(l-oos9)   -  Pr[l-coa(9-<x)]  ]  oos9d9  -  0. 
Collecting  terms,   we  get 

"EI"     a  [Jo^"rSi"  9  +  Kroos9  "  Kroos29)d9  +  JK (rsin29  +  Kroos9  -  Kroos29)de] 

+  "EI"     a  [/    (-sinec°s9  +  Ksin9cos9)d9  +  j     (-sin9oos9  +  Ksin9co89)d9] 

+  -EJ— Ta  [J"(sin29  +  Kcos2e)d9  +J<x(sin29  +  Koos29)d9] 


"EI""  J       [frain(9-«)sin9  -  Prcos9  +  Prcos(9-a)cos9]d9. 


-(42) 


E1-  (40),  (41)  and  (42)  can  be  written  in  a  matrix  form  as 

„Ft 


(Q)  \vF 


G 


.  a 


-(43) 


T 
where  G,  H,  S   corresponds  to  the  vector  D  in  Eq.  (35). 

Comparing  Eq.  (43)  with  Eq.  (35),  we  know  that 

(Q)  ■  (A) 


-(44) 


So  we  get 


r  Fi 
F 
a 

M1" 
a 

_P 

T 

L  a 


-  -(A)"1 


0 

H 

LsJ 


-(45) 
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where 


G  **  "EI""  f   [-prsin9sin(9-c«)  -  KPr 
x       °< 

+  KPrcos(e-a)  +  KPrcos9  -  KPrcos(9-«)cos9]d9, 

H  =  "EI"  I   [-Prsin(e-«)cos9  -  KPrsin9  +  KPrco3(9-«)sinGld9, 

S  =  "EI™  J   [Prsin(9-a)sin9  -  Proos9  +  Prcos(9-«)cos9]< 
Integrate  and  get 


d9. 


5 


Pr      t 
G  "  ~4EI]   (1-K)°OS0((sin2^-3in2a)   +  (K-l)sina(cos8f-cos2«) 

-   2(4'- «)     [(K+l)oos«+  2k]      +  4Ksin(£-«)  +  4K(sin£-sin<x)], 


-(46) 


H  "  -2EI"[2^"0(^K+1)3ino<+  4K(cos^-oos«)   -  co8«(cos2f-cos2«)(K-l) 
r 

-sin«(3in2^-sin2«)(K-1 )]    ,   (4y) 

S  "  "ZEI"   [2(4'-0<)(K+1)cos«-  4K(8in$-sin«)   +  cos«(sin2£-3in2«)(K-l) 
r 


-8in«(cos2f-cos2o<)(K-l)] 
Rewrite  Ec[.(38)   and  get 


-(48) 


l9. 


rsin9 


0  0 

oos9       -sin9 


\r(l-oos9)       sin9  oos9/         T 


r   P 
rF 


.  MF  •    - 


Prsin(9-«) 
Pr[l-cos(9-«)] 


-(49) 


The  fixed-end  stresses  at  end   "b"  can  be  obtained  by  substituting  the 
values  of  If,   F,   r   into  Eq.(49),   and  changing  9  to  j>.       Thus  we  get 


„Fn 


K 


L   bJ 


0  0 

cos<£       -sin^> 

^r(l-cos£)       sin^>        cos-£/ 


(A) 


-1 


'  G 
H 

LsJ 


+  <j  Prsin(<£-oc) 
Pr[l-cos(<j>-oOl 


-(50) 
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STIFFNESS  MATRIX  OF  A  STHAIGIIT  HEMBEfi 

A  straight  beam  with  end  "a"  free  and  end  "b"  fixed  is  shown  in 
Fig.  7.  Assume  that  the  loading  and  deformation  are  all  in  the  Z  direction. 


.; 
-  ■ 

1 

a 

Ml 
f 

n; 

' 

A 

T| 

$ 

i 
) 

| 

h 

X 

L 

Fig.  7.  Straight  beam  with  one  fixed  and  one  free  end 

The  properties  of  segment  ab  can  be  described  by  the  partitioned  matrix 
equation, 


$ab 
£ga 


bO 


-(51) 


where 


-FF] 
a 

>\  W  ■ 
a 

p 

r 

a 

[■1*1 

a 

•  M*"  • 

a 


and 


ya 


zb 

yt> 

LRxb- 
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Dzi  A.  and  Rx  indicate  the  displacement  in  the  Z  direction  and  rotations 
about  Y  and  X  axis  respectively. 

Then  them's  in  the  stiffness  matrix  can  be  obtained  as  follows. 
At  any  section  of  a  distance  x  from  end  "a"  the  stresses  can  be  expressed 

as 


i* 


T 
u  x 


1        0       0 

rFF] 

a 

x       1        0 

a 

0       0        1/ 

F 
T 
u    aJ 

-(52) 


The  total  strain  energy  stored  in  the  beam  due  to  the  acting  forces, 
neglecting  the  effect  of  vertical  shear,  will  be: 


Jo 


H2dx 


2 

,,   T  dx 

+  /    X 
'  c 


>o     2EI   T  )  0   ~SGJ~~  ' 

-  A  [^  -  «W  ♦  -K>V  ♦  KL(/)2]. (53) 

By  Castigliano's  theorem 


°*a  "  5  -  "A"  W  *  -|-3<]  • 


ya 


"5  =  "A-  K  -  *fl  •  - 


-8U 


"SET  2KL< 


-(54) 

-(55) 
-(56) 


E1-(54),  (55)  and  (56)  can  be  expressed  in  the  form  of  a  matrix  equation 
r.3        ,2 


ya 


L 

"sr 


r  F 
a 

a 


-(57) 
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r  F " 
F 
a 

/    12EI 

6EI 

-_-r 

•  \ 

<- 

a 

>        = 

6EI 

4EI 
L 

0 1 

P 

\      o 

0 

GJ    / 

"VI 

So  we  get 

/    12EI 

6EI 

0       \ 

#ab     = 

6EI 

4EI 
"E" 

0 

I     o 

0 

GJ    / 
"VI 

ya 


-(58) 


-(59) 


To  get  the  stresses  at  end  "b",  substitute  the  value  of  Eq.(58)  Into 
Eq.(52)  and  let  x-L.  Then 


[F" 

/1 

0 

0 

"d 

za 

< 

kJ 

>.  —  _ 

-I, 

o 

1 

0 

0 

l] 

4>h>     < 

R 
ya 

B 

-   xa 

-(60) 


So  we  get 


V 


k 


12EI 


j>%a 


6  EI 


L* 

_..r 

0 

6EI 

2EI 

-.r- 

-  — r- 

.     0 

0 

0 

GJ 

"I 

12EI 

est 

l' 

~?- 

0 

6EI 

2EI 

IT 

-  ..L- 

0 

0 

0 

GJ 
1" 

<  B 


.ya 


-(61) 
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To  determine  the  values  of  ^>Ba  and  4'^tl>  we  keep  the  end  "a"  of  the 
beam  ab  in  Fig.  7.  fixed  and  let  the  end  "b"  free. 

At  any  section  of  a  distance  x  from  end  "b"  the  stresses  can  be 
expressed  as  in  Eq.(62),  or 


T 
U  x 


<  +  f£*. 

V 


-(62) 


T 
.  x 


o\ 

r    F  1 

0 

X 

i/ 

F 

KJ 

-(63) 


Similarly  following  the  steps  from  Eq.(53)  to  Eq.(57),   we  can  get 


ssb 


<  B 


yb 


xbJ 


5EI" 


L 

~2El~ 
L 


0 
0 

"sr/ 


< 


"b-l 


-(64) 


K 

i 

< 

k\ 

\ 

12EI 
..... 

6EI 
...r 

0 


6EI 

...r 

4EI 
0 


GJ 

•"V 


zb 


yt 


xbJ 


-(65) 


Thus 


£fea 


12EI 
....3 

6EI 

-  _..r 

0 

6EI 
"2" 

4EI 

0 

0 

0 

GJ 

-(66) 
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To  get  the  stresses  at  end  "a",  substitute  the  value  of  Bq.(65)  Into 
Eq.(6j)  and  let  x=L,  then 

>F1  /  1         0         0 


\o 


L  1  0 

0  ,/ 


zb 


£k  \  v  ■ . 


.   xb 


■(67) 


Thus 


/12EI 

6EI 

hr 

.  _.-r 

0 

n*  .- 

6EI 

___r 

_2EI_ 

0 

\      o 

0 

GJ 

"T 

-(68) 
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NUMERICAL  EXAMPLE 

The  gridded  framework,  with  straight  members  and  members  curved  in 
plan,  shown  in  Fig.  8.  is  fixed  at  points  1,4,5  and  8.  Points  2, J, 6  and  7 
are  free  to  rotate  and  to  translate  perpendicular  to  the  plane  of  the  frame. 
The  frame  could  be  girders  of  a  bridge.  The  cross  sections  A-A,  B-B  and  C-C 
of  Fig.  8.  are  shown  in  Fig.  9  and  Fig.  10  respectively.  The  analysis  of 
stresses  will  be  separated  into  two  parts.  One  is  the  calculation  of 
stresses  due  to  dead  load  and  the  other  the  calculation  of  stresses  due  to 
live  load.  When  we  do  the  calculation  of  stresses  due  to  live  load,  the 
method  of  influence  lines  will  be  used  in  this  report. 


Fig.  8.  Plan  view  of  a  frame  with  curved  and  straight 
members 
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Fig.  9.   Section  A-A  in  Pig. 


1-4" 


Fig.  10.   Section  b-B  and  C-C  in  Fig.  8. 
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The  following  notation  and  assumptions  are  used  in  the  calculations: 


Modulus  of  elasticity  of  steel 
Ratio  of  E  to  E 


Bb  -  29,000,000  lb/in 

n  .  9 

E  .  3,220,000  lb/in2 


Modulus  of  elasticity  of  concrete 
Moment  of  inertia  of  curved  beams 
Moment  of  inertia  of  straight  beams 
'niform  desd  load  for  curved  members 
Uniform  dead  load  for  straight  members   q  =1.24  kips/ft 
The  torsion  factor  of  curved  beams  is* 


Ir  -  863,370  in4 
I.  -  228,750  in4 
qo  -  1.70  kips/ft 


Jo fj?       (L2)5[54)l   _  .  1J7>000  in4. 

0   3.58(b2  +  d2)    3.58[(22)5  +  (54)5] 

The  torsion  factor  of  straight  beams  is 


If 


bV 


»3, 


„ S- (I6)l46)f._  _ 

5.58(b2  +  d2)    3.58[(16)2  +  (46)2] 


2 — 72T  "  ~Z'ZZ7t"X^'~" — 2i"  "  47'000  in 


-4-   -  2.25*. 

where     E  «  modulus  of  elasticity  in  moment, 

G  »  modulus  of  elasticity  in  shear, 
then  the  value  of  K  for  curved  beams  is 


EI 

r 

"5J- 


2  2S  863370 
*,Z5  "37505" 


14.2. 


Refer  to  reference  11  in  the  Bibliography. 


DETERMINATION  OF  STRESSES  DUE  TO  DEAD  LOAD 

The  Slope-Deflection  procedure  in  the  computer  method  of  stress  analysis 
will  be  used  in  this  report.  First  the  fixed-end  stresses  will  be  determined 
in  part  A.  Then  these  fixed-end  stresses  will  be  used  to  evaluate  the 
member  end  stresses  in  part  B. 

A.  The  calculation  of  the  fixed-end  stresses: 

In  this  example  we  have  y»  -~-   and  K  ■  L4.2  for  curved  members. 

To  evaluate  the  fixed-end  moments  and  torques  in  Eq.(32)  and  (33)  for 
each  curved  member,  we  shall  use  the  electronic  computer.  FORGO  source 
program  1  is  used  for  this  purpose  and  is  shown  on  pages  61  and  62.  Answers 
appear  below  on  this  page. 

For  the  convenience  of  computer  calculation,  some  substitutions  have 
been  made  in  the  source  program.  The  values  BMF  and  TMF  in  the  source 
program  are  used  instead  of  the  values  Ma  and  Tft  in  Eq.(32)  and  (33). 
A,B,C,D,  and  E  are  the  constants  in  Eq.(29).  We  shall  use  q  as  a  positive 
value  since  we  have  considered  that  q  is  negative  as  we  derived  Eq. (32) 
and  Eq.(33). 

ANSWERS  OF  PROGRAM  1 

A  It  C  D  E 

5. 3C58C110   1.5441067Ci   4.?4?3<J8?0    .^1113570   -.01567320 


PfffiING-'""'Hll  TCRSICN-f-'CN'ENT  P 

-75.369'  '.  5<    i                            -.34668919  63.661976C0 

'                                B                               C                                D  F 

5.3>./58C110       1.54410670  4.24239830          .01103570  -.01567320 


8END1NG-M0MENT  TCRSICN-KCMENT  R 

-  59.7U57240(.  -.27463982  56.6617/600 
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Now  with  the  data  given  in  the  answers  of  program  1 ,  the  fixed-end 
stresses  for  each  member  can  be  listed. 

In  the  list  we  use  the  notation  m£  which  will  indicate  the  fixed-end 
moment  of  member  (1)  at  end  2(the  first  subscript  refers  to  the  end  of  a 
member  and  the  second  subscript  refers  to  the  number  of  bars)  and  the  same 
subscription  will  be  used  for  fixed-end  torques  and  forces. 

Because  members  (l),  (2)  and  (3)  have  the  same  curvature  and  length, 
the  fixed-end  stresses  due  to  uniform  load  should  be  the  same.  So  we  get 


M^1   =  -  75.369  ft-kips, 

*- 

+  75.369  ft-kips, 

T^1  =  -  0.347  ft-kips, 

_F 

T 

-  0.347  ft-kips, 

F^  -  +  19.778  kips, 

*♦ 

+  19.778  kips; 

"22  =  "  75*369  ft-kips, 

4- 

+  75.369  ft-kips, 

TF22  =  -  O.347  ft-kips, 

_F 

T       =. 
32 

-  O.347  ft-kips, 

P^2  -  +  19.778  kips, 

F 
F        ■ 
32 

+  19.778  kips; 

M^j  -  -  75.369  ft-kips, 

43 

+  75.369  ft-kips, 

T^     -  -  O.347  ft-kips, 

_p 

T       . 
43 

-  O.347  ft-kips, 

F33  -  +  19.778  kips, 

y 
P       = 
43 

+  19.778  kips. 

Because  members  (6),  (7)  and  (8)  have  the  same  curvature  and  length, 
the  fixed-end  stresses  due  to  uniform  load  should  be  the  same.  So  we  get 

+  59.706  ft-kips, 
-  O.275  ft-kips, 
+  17.603  kips; 


«i- 

-  59.706  ft-kips, 

<e 

-  0.275  ft-kips, 

4 

^6  = 

+  17.603  kips, 

4 

J3 


MC7  -  -  59.706  ft-kipo,  "77  "  +  59.706  ft-kips, 

Tgy  -  -  O.275  ft-kipB,  T^?  .  -  0.275  ft-kips, 

Fg?  -  +  17.603  kips,  F?  .  +  17.603  kips; 

1^8  -  -  59.706  ft-kips,  M^8  .  +  59.706  ft-kips, 

T^Q  =  -  0.275  ft-kips,  TgQ  -  -  0.275  ft-kips, 

F?8  -  +  17.603  kips,  Pg8  .  +  17.603  kips. 

Because  members  (4)  and  (5)  have  the  same  length,  the  fixed-end  stresses 
should  he  the  same.  So  we  get 

Mj4  -  +  5.06  ft-kips,  M^  +  -  5.O6  ft-kips, 

Pg4  -  +  4.34  kips,  F^4  =  +  4.34  kips; 

M^5  =  +  5.06  ft-kips,  mL   -  -  5.O6  ft-kips, 


F 
EL  -  0, 


^5  ■  °'  -75 

F^5  -  +  4.34  kips,  ^75  "  +  4,J4  kiP8' 


J4 


B.  Use  the  determined  fixed-end  stresses  to  evaluate  the  member  end  stresses. 

This  part  of  the  calculation  will  be  separated  into  6  steps. 
1.  General  approach. 

Define  the  statics  matrix  A  expressing  the  balancing  forces  {q}  at 
the  joint  in  terms  of  the  member  end  stresses  {p}  .  Then  by  definition, 

{ft]  -  [A]  {P}  .  (69) 

Define  the  geometry  matrix  B  expressing  values  of  elastic  member  end 
rotations  and  displacements  {e}  in  terms  of  joint  rotations  and  displace- 
ments {X}  .  Then  by  definition, 


where 


{e}  =  [Bj  {X}  ,   (70) 

[B]   "  [A]  T (71) 

Define  the  stiffness  matrix  K  expressing  the  member  end  stresses  Ip} 
in  terms  of  the  elastic  member  end  rotations  and  displacements  {e}  . 
Then  by  definition, 

{?}  -  [K]  {e}   .  (72) 

Substituting  Eq.(70)  and   (71)   into  Eq.(72), 

{P}      =      [K]    [A]T{X]     .      [G]{X) (73) 

Substituting  Eq.(73)   into  Eq.(69), 

JQj      =■      [A][K][A]T{X]      .      [H]{X} (74) 

Therefore 

lx}    "     [Hi"1  W ,- (75) 
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2.  Draw  the  Q-X  and  P-E  diagrams. 

In  Fig.  11,  the  Q-X  diagram  shows  the  joint  balancing  forces  and 
corresponding  deformations.  All  the  forces  are  shown  in  the  positive 
direction. 

In  Pig.  12,  the  P-E  diagram  shows  the  member  end  stresses  and  corres- 
ponding member  end  deformations.  To  keep  the  picture  clear  only  the  end 
stresses  and  corresponding  deformations  at  end  2  of  members  (l),  (2)  and 
(4)  are  presented.  Actually  at  each  end  of  each  member  there  will  be  three 
pairs  of  stress-drformation  expressions.  All  the  member  end  stresses  are 
3hown  in  the  positive  direction.  Also,  since  members  are  perpendicular  to 
each  other  at  joints  the  torque  of  the  straight  bar  is  in  the  direction  of 
the  moment  of  the  curved  bar. 

In  Fig.  13,  the  diagram  shows  the  equilibrium  condition  of  balancing 
forces  and  end  stresses  at  joint.  According  to  this  condition  the  sum  of 
balancing  forces  and  end  stresses  at  a  joint  has  to  be  equal  to  zero. 
Therefore,  the  balancing  forceB  will  be  the  sum  of  the  end  stresses  at 
member  ends.  For  example, 

F2  =  F2  +  F2  +  F2'  (?6) 

M2  -  k\   +  l\   +  M*,  (77) 


at  joint  2  in  Fig.  13. 

It  is  noted  that  the  values  to  which  F„,  M„,  .  .  .  ,  refer  in  Eq.(76) 
and  (77)  are  the  end  stresses  at  the  joints.  They  are  not  the  member  end 
stresses.  The  drawing  of  the  equilibrium  diagram  is  an  important  procedure 
by  which  we  can  establish  matrix  A  in  Eq.(69). 
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Pig.    11.      Q-X  diagram 


Pig.    12.      P-E  diagram 


57 


joint  2 


joint  3 


joint  6 


joint  7 


Fig.  13.  Equilibrium  diagrams  of  joints  2,3,6  and  7 
3.  Establish  the  matrices  A,B  and  K. 

To  establish  matrix  A,  we  inspect  every  joint  and  write  the  equilibrium 
equations  as  in  Eq.(76)  and  (77)  for  each  balancing  force  component.  Then 
matrix  A  can  be  written.  For  example,  the  entries  of  the  first  row  in  Eq. 
(78)  are  all  zero,  except  the  fourth,  seventh  and  tenth  which  are  equal  to 
one.  So  we  can  write  Fg  ■  T?   +  FJ  +  K.  This  is  also  what  Eq.(76)  indicates. 
Matrix  B  is  simply  the  transpose  of  matrix  A. 

If  we  partition  the  matrices  A  and  B,  we  can  rewrite  them  in  the  following 
form 


[A]   =  [A1  |A2], 


[Bj 


Bl 
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H: 


z2 

R2 
Rr2 

"t2 
r3 


E' 


t3 

z3 


r3 


t3 
z3 

z4 

& 

r4 

»3 


\l  \l  Rt2Dz2  Rr2  Rt2Dz2  Rr2Rt2  Dz2Rr2  Rt2 
"o.    0.    0.    0.    C.    0.    0.    0.    0.    0.    0.    o«" 

v.   o.   o.  o.   o.   o.   o.   o.   o.  n.  o.   o. 

n.     0.    0.    0.    0.     0.    0.    0.    0.    0.    0.    0. 

i .  c.  n.  o.  o.  o .   o ,  o.  o.  o.  o.  o. 

A.  l.  0.  n.  o.  o.  o.  o.  o.  o.  o.  o. 

«.  n..l.  o.   n.  o.  o.  o.  o.  o.  n.  o. 

i.  o.  o.  n.  o.  o.  o.  o.  o.  o.  o.  o. 

r\,  o.    ).   n.   o.  o.   o.   o.   o.   o.  o.   c. 

".  l.   o.  o.   o.  o.   o.  o.  o.  o.  o.  o. 

1.    0.    0.    n.    0.    0.    0.    0.    0.    0.    0.    0. 

n.    i.   o.   n.   o.   o.   o.   o.   o.  o.  o.   n. 
o.   o.   ].   n.   o.   o.   o.  o.   o.   o.  o.   p. 

0.  0.  0.  1.  0.    0.  0.  0.  0.  0.  0.  0. 

n.  0.  0.  0.  1.    0.  0.  0.  0.  0.  0.  0. 

^.  0.  0.  0.  0.     1.  0.  0.  0.  0.  0.  0. 

n.  0.  0.  ].  0.    C.  C.  0.  0.  0.  0.  0. 

-i.   o.   o.  n,  o.   i.  o.  o.  o.  n.  o.  o. 

".     ".    0.    o.     1.     C.    0.    0.    0.    0.    0.    0. 
".     0,     ".     1.    0.     0.    0.    0.    0.    0.    0.    0. 

".    ".   o.   o.    l.    o.   n.   o.   o.   n.   o.   o. 

".    0.    0.    0.    0.     1.    0.    0.    0.    0.    0.    0. 
".    0.    0.    0.    0.    0.    0.    0.    0.    0.    0.    0. 

".   o.   o.   n.  o.   o.  o.  o.  o.  o.  o.   o. 


r'     o.  o.  o.  n.  o.  o.  o.  o.  o.  o.  o.  o. 

t4  ^  -1 


-(80) 
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M 


°%6 


z6 


r7 

Kt7 

z7 

«\ 
i'7 

B5 

Rt7 
0 
z7 
8 
r7 
8 
t7 


z8 


i-8 


Dz2Rr2  Rt2Dz2  Rr2  Rt2Dz2  Rr2Rt2  Dz2Rr2  Rt2 

"n.  0,  r.  o.  o.  c.  r.  o.  o.  o.  o.  o. 

".  o.  o.  o.  o.  c.  c .  o.  o.  o.  o.  o. 

0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0. 

n.  0.  0.  0.  0.  0.  1.  0.  0.  0.  0.  0. 

a.  n.  o.  o.  o.  o.  o.  i.  o.  o.  o.  o. 

r.  0.  0.  0.  0.  0.  0.  0.  1.  0.  0.  0. 
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"t8U 
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-(81) 
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A.,  A2,  B.  and  B,  all  represent  the  matrices  and  they  are  shown  in 
Eq.(78),  (79),  (80)  and  (81 )  respectively.   It  is  noted  that  the  notation 
of  forces  or  deformations  listed  on  the  top  and  on  the  left  side  of  the  matrix 
is  not  included  in  the  matrix.   It  is  only  for  reference.  For  example, 
matrix  A  is  a  member  end  stress  and  balancing  forces  transformation  matrix. 
So  we  list  the  notation  of  member  end  stresses  on  the  top  and  balancing 
forces  on  the  left  side  of  the  matrix  to  show  to  which  pair  of  stresses  and 
forces  each  entry  of  the  matrix  refers. 

To  establish  the  matrix  K,  we  first  have  to  determine  the  submatrices 

A  A 

K.,  K„.,  .   .   .,  etc.  in  Eq.(24). 


The  values  of  K   and  K   are  determined  in  program  2  which  starts  on 

2    -3 

22  "  A3J' 


page  62.  Answers  appear  on  page  43.   In  this  example  K^  »  K£„ 


Md  K55  "  K66  "  K77- 

The  values  of  K   and  Kg  are  determined  in  program  3  which  starts  on 

page  65.  Answers  appear  on  page  44.   In  this  example  K._  -  K.  -  k'  , 

and  K656   -  K}7   =  K9le 

The  values  of  K  _  and  K,g  are  determined  in  program  4  which  starts  on 


page  67.  Answers  appear  on  page  44.   In  this  example  K„_  -  K 


2  -«? 


22   "33   "44' 
-d  4  "  K77  ■  K88- 

The  values  of  K   and  K,^  are  determined'  in  program  5  which  starts  on 

page  69.  Answers  appear  on  page  44-   In  this  example  K„.  •  K,„  -  K?.. 

and  K665  .  K]6  .  K8er 

The  values  of  K^2  and  Kjjg  are  determined  in  program  6  which  starts  on 

page  71.  Answers  appear  on  page  45.   In  this  example  C,  -  K^   and 

K4  .  K5 
K26   K37" 

The  values  of  Kgg  and  icL  are  determined  in  program  7  which  starts  on 
page  73.  Answers  appear  on  page  45.  In  this  example  IQ,  -  SL,  and 


43 


K4     .  K5 
K62  "  K73' 

After  getting  these  submatrices,   we  can  insert  them  into  Eq.(24)   and 
get  the  matrix  K. 

For  convenience  in  writing  the  source  program,   some  changes  have  been 
made  in  the  notations  in  the  derived  formulas  as  they  are  used   in  the 
source  programs. 

In  program  2,j>  ,   K  and   [a]"     in  Eq.(9)  have  been  changed   into  X,   GK 
and   [Hj ,    respectively. 

In  program  3,  <$>  and  K  in  Eq.(l2)  have  been  changed   into  X  and  OK, 
respectively. 

In  program  4,$,   K,    [b]   and   [b]"1    in  Eq.(l6)  have  been  changed   into 
X,   GK,    [A]    and   [h] ,   respectively. 

In  program  5,4>,   K  and    [b]"     in  Eq.(l8)  have  been  changed   into  X,   GK 
and  [A]      ,   respectively. 

The  answers  of  program  2-7  are  listed  as  follows: 


AMSWFRS    CF    PROGPAM    2 


STIFFNFSS-MATRIX     INV.(A)     IN    EO.I9) 


2C44Q.9P00 

■22705'.. 2800 

1 t 107.4900 


-227054.2800 
1266355.0000 
-278175.0500 


14197.4900 

-27«175.0500 

92871.4480 


63.661976 


11 


STTFFNF^-MATRIX     INV.(A)     IN    FO  .  (  9  ) 


2900?. 9020  -286620.0200  17922.0810" 

■286620.0300  3669879.6000  -312540.7200 

17922.0810  -312540.7200  104344.7600 


56.661976 


55 


ANSWERS  OF  PROGRAM  3 

STIFFNFSS-MATRIX  IN  FQ.I12) 


44 


-20449.9000      227054.2800  -14197.4900 

-227052.5100     1779295.9000  -15968.6430 

-14197.2000       15966.6400  -46637.3990 

STIFFNESS-MATRIX  IN-EQ.U2) 


63.661976 


12 


t  -2°003.9020 
286617.7800 
-17021.7130 


286620.0200 

1999109.6000 

17939.1100 


56.661976 


-17922.0810 
-17941.4300 
-52398.9710 


"K56 


ANSWER5  OF  PROGRAM  4 

STIFFNFSS-MATRIX  IN  EQ.I16) 


20449.9000      227054.2800       14197.4900 

227054.2800     3266355.0000      278175.0500 

141  )7.4900      278175.0500       92871.4480 


63.661976 


■Kl22 


STIFFNESS-MATRIX  IN  F0.I16) 


29013.9020      286620.0200       17922.0810 

286620. C300     3669879.6000      312540.7200 

17922. C810      312540.7200      104344.7600 


56.661976 
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ANSWERS  3F  PROGRAM  5 

STIFFNFSS-MATRIX  IN  EQ.U8) 


-2^449,9000 
227052.5100 
-14197.2000 


-227054.2800 

1779295.9000 

-15966.6400 


-14197. 4900v 

15968.6430 

-46637.3990 


63.661976 


K21 


STIFFNFSS-MATRIX  IN  F0.(18) 


-29003.9020 
286617.7800 
-17921.7130 


-286620.0200 

1999109.6000 

-17939.1100 


-17922.0810^1 

17941.4300 
-52398.9710 


,K6 


56.661976 


ANSWERS  OF  PROGRAM  6 


STIFFNESS-MATRIX  FCR  NEAR  FNO 

178954.0800      626339.7900  0.0000' 

676339.7900     2922916.7000  0.0000 

o.onno           0.0000  66728.3960 
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STIFFNESS-MATRIX  FCR  FAR  END 

-178954.0800     -626339. ?900  -0.0000' 

626339.7900     1461458.^000  -0.0000 

-0.0000                          -0.0000  -66728.3960. 
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ANSWERS    OF    PROGRAM    7 

STJFFNFSS-MATRIX    FOR    NEAR    FNO 
178954.0800  -676339.7900  0.0000" 

-6763:9.7900  7977916.7000  0.0000 

0.0000  0,0000  66728.3960. 


STIFFNESS-MATRIX    FOR    FAR  FND 

-178954.0800              626339.2900  -0.0000" 

-6763T9.7900            1 461458 . 3000  -0.0000 

-0.0000                          -0.0000  -66728.3960 


K4 


K4 
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4.      Calculate  matrices  G  and  H~ 


Through  matrix  multiplication  and  inversion,  the  matrices  G  and  H~ 
can  be  easily  evaluated.  Programs  8  and  9  shown  on  pages  74  and  75  are 
examples  of  matrix  multiplication  and  inversion. 


5.     Determine  the  vector 


Using  the  fixed-end  stresses  at  member  ends  liBted  on  pages   32  and 
33.   vector  Q  can  be  determined  from  information  shown  in  Fig.    14. 
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■  19.778 


+75.360 


■19.778  ^-r.4,34  "^5*369^   ,Q 


+0.347 

+0.347^^ 

Tl — Oy_5.06 


K2  =  -43.806  kips     M2  =  0 


T2  =  -4.366  ft-klps 


-10.775 


+75.360 

o 


+0.347 


--.-^<r-4.34     -t^5<?10       -^fC^-5. 


06 


F„  =  -43.806  kips      M„  =  0 


T„  =  -4.366  ft-kips 


-4.3  4 


-17.603 


17.603 


.ik     +50.706 


6  -59. 706^0^° 


+  5.06 
"Z^      +0.275 

+  0.275>^v— — ~) 

6     ' 


Ff)  =  -30.546  kips      MR  =  0 


Tfi  =  +5.61  ft-kips 


-4.34 


-17.603 


17.603 


>>.     +59.706 


-59.706 


r> 


+5.06 

^     +0.275 

+0.27^T     ' 


F?  =  -39.546  kips      M?  =0  T?  =  +5.61  ft-kips 

Fig.  14.   Balancing  forces  at  Joints  2,3,6  and  7 


Therefore  the  vector  Q  is 
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W 


'V 

-43.896 

M2 

0.0 

T2 

-  4.366 

F5 

-43.896 

M3 

0.0 

F6 

•  .  < 

-  4.366 

-39.546 

\ 

0.0 

T6 

+  5.61 

F7 

-39.546 

«7 

0.0 

-T7- 

.+  5.61  . 

6.  Calculate  the  member  end  stresses  as  indicated  in  vector  P  in  Eq.(73). 

Substitute  the  value  of  matrix  H~  and  vector  Q  into  Eq.(75).  Vector 
X  can  then  be  determined.  Then  substitute  the  value  of  matrix  G  and 
vector  X  into  Eq.(73).  The  vector  P  can  be  thus  determined. 

Here  we  have  used  a  source  program  of  matrix  times  vector.   It  is 
shown  in  program  10  on  page  76. 


Vector  P  was  partitioned  in  the  form  of   (PI 


Combining  the  fixed-end  stresses  with  the  member  end  stresses,  vector 
P,  we  get  the  final  member  end  stresses  which  are  listed  in  Tables  1  and  2 
on  pages  85  and  86. 
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DRAW  INFLUENCE  LINE 

In  this  section,  we  shall  draw  the  influence  line  for  each  stress 
component  of  each  section  shown  in  Pig.  15.  If  more  seotions  were  to  be 
taken  the  influence  lines  would  be  smoother.  However,  we  will  define  the 
ordinates  only  at  sections  shown  in  Fig.  15  when  we  draw  the  influence 
lines  later.  The  change  of  values  between  seotions  will  be  assumed 
linear. 


Fig.  15.  Plan  view  of  a  framework  with  curved  and  straight 
members 

First  the  fixed-end  stresses  will  be  determined  in  part  A.  Secondly 
the  influence  line  will  be  drawn  in  part  B. 

A.  The  calculation  of  the  fixed-end  stresses. 

Before  determining  the  fixed-end  stresses,  the  constants  Q,H  and  S  in 
Eq-U6),  (47)  and  (48)  were  evaluated  in  program  11  starting  on  page  77. 
Answers  appear  on  page  49.  Then  substituting  the  values  of  Q,  H  and  S  into 
Eq.(45),  we  get  the  fixed-end  stresses  at  end  "a".  They  are  determined  by 
program  12  starting  on  page  80.  Similarly  if  the  values  of  Q,H  and  S  are 
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substituted  into  Eq.(50)    the  fixed-end  stresses  at  end   "b"  can  be 
determined.      This  is  performed  by  use  of  program  13  starting  on 
page  82.      Answers  of  program  12  and  13  appear  on  page  50  and  51. 


AmSVFR?  CF    PROGRAM  ]1 

R     63.66197600 


-.000140844240   -.000009378848   -.000006671969 
R     63.66197600 


Z 
.08726646 


G  H  S  Z 

.17453293 

-.000066033581   -.000003845325   -.000001600150 
R     63.66197600 


-.000017110796   -.000000864254   -.000000061282 
R     56.66197600 


-.000099305638   -.000007429725   -.000005285392 
R     56.66197600 


-.000046558573   -.000003046185   -.000001267605 
R     56.66197600 


Z 
.26179939 


Z 
.08726646 


Z 

.17453293 


-.000012064381   -.000000684644   -.000000048546 

STOP   END  CF  PROGRAM  AT   STATEMENT  0060  +  00  LINES 


Z 
.26179939 
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ANSWERS  CF  PROGRAM  1? 


63.66197600 
FIXFD-FND  STPFSSFS 


.84546821 

-3.20061630      .01030862 

.08726646 

R 

63.66197600 

FIXFD-FND  STRESSES 

F 

.50000071 

R 

M              T 
-2.87813180      .01644585 

63.66197600 

Z 
.17453293 

FIXFD-FND  STRESSES 


.15455155 

-1.07916620      .00820781 

.26179939 

R 

56.66197600 
FIXED-END  STRESSES 

F 
,8454683? 

M              T 
-2.84868920      .00917506 

Z 
.08726646 

R 

56.66197600 

.50000081 
R 


FIXFO-FNn  STRPSSPS 


li 

•2.56166610 


.01463761 


56.66197600 
FIXED-END  STRESSES 


.17453293 


F  M  T  Z 

.15455148     -.96050480      .00730520      .26179939 
STOP   FND  CF  PROGRAM  AT   STATEMENT  0050  +  00  LINES 


ANSWERS  CF  PR^GPAM  13 


61.66197600 

FTxr:n-FNn  strcssfs 
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.1  ^"i^lflO 
R 


M 
1  .07909800 


.00822890 


63.66197600 
FIXED-FND  STRESSES 


.08726646 


.49999930 

2.87811500      .01644997 

.17453293 

R 

63.66197600 
FIXFD-FNO  STRESSES 

F 
.84544850 

M              T 
3.20024310      .01027509 

Z 
.26179939 

R 

56.66197600 

FIXED-END  STRESSES 


.15453170 

.96044500      .00732360 

.08726646 

R 

56.66197600 

FTXFD-FND  STRFSSFS 

F 
.4999  >920 

M              T 
2.56164840      .01464101 

Z 
.17453293 

R 

56.66197600 

FIXED-END  STRESSES 

F  M  T  Z 

.84544860     2.84835840      .00914525      .26179939 
STOP   END  OF  PROGRAM  AT   STATEMENT  0060  +  00  LINES 
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B.  Draw  the  influence  line. 

This  will  be  separated  into  four  steps. 

1  ■  Determine  the  vector  Q. 

Using  the  fixed-end  stresses  determined  in  programs  12  and  13,  we  can 
follow  the  same  procedure  to  determine  the  vector  Q  under  a  unit  load  as  we 
determined  the  vector  Q  under  uniform  load  on  page  46.  For  every  acting 
position  of  the  unit  load,  there  will  be  one  corresponding  vector  Q.  Results 
are  listed  in  Table  3-8  starting  on  page  87. 

2.  Calculate  final  member  end  stresses. 

Substituting  the  vector  Q  into  Eq.(75),  we  get  values  of  vector  X. 
Then  from  Eq.(73),  we  get  values  of  vector  P,  the  member  end  stresses. 
Combining  the  member  end  stresses  with  the  fixed-end  stresses,  we  get  the 
final  member  end  stresses.  They  are  listed  in  Table  9-12  starting  on  page 
90. 

3.  Calculate  the  stresses  at  a  section  between  joints. 

The  stresses  at  a  section  between  joints  can  be  determined  by  substi- 
tuting final  member  end  stresses  into  Eq.(l0),(l3)  and  (52).  Results  are 
listed  in  Tables  13»14>15  and  16,  in  which  the  notation  referring  to  each 
stress  component  is  as  it  is  defined  earlier  in  the  report.  For  example, 
F18  simply  represents  the  vertical  force  at  section  18  shown  in  Fig.  15. 
However  here  we  did  not  use  any  number  of  bar  as  the  subscript  with  the 
notation  as  in  Table  13-16  starting  on  page  98. 

4.  Draw  the  influence  line. 

When  we  draw  the  influence  line  we  indicate  only  the  ordinate  above 
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the  beam.  Although  data  are  available  to  draw  an  influence  for  each  stress 
component  of  each  section,  we  will  draw  only  the  influence  line  for 
section  13  as  an  example.  For  the  convenience  of  drawing  the  influence 
line  all  the  curved  beams  have  been  lengthened  and  have  become  imaginary 
straight  beams. 

Figures  16,  17  and  18  show  the  influence  lines  of  vertical  force, 
moment  and  torsion  at  section  13,  respectively.  The  values  of  coordinates 
of  these  lines  cp.n  be  found  from  Tables  13,  14,  15  sjid  16.  For  example, 
to  find  coordinates  of  the  influence  line  for  vertical  force  at  section  13, 
we  turn  to  page  99  where  is  the  Table  1J.   In  the  left  most  column  of  this 
table,  we  find  F.  .  Then  values  appearing  in  the  same  row  with  F...  are 
vertical  forces  at  section  1 3  due  to  a  unit  load  acting  at  each  different 
section  in  turn. 
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Fig.  16.   Influence  line  for  vertical  force  at  section  13 
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Fig.  17.   Influence  line  for  moment  at  section  13 
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Kig.  18.   Influence  line  for  torsion  at  section  13 
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CONCLUSION 

All  the  formulas  derived  in  this  report  are  intended  for  general 
application  to  any  kind  of  structure  with  members  which  are  straight 
or  curved  in  plan.  Any  more  complicated  structure  can  be  analysed  by 
the  same  procedure.  In  case  the  number  of  unknowns  in  the  stiffness 
matrix  exceeds  the  capacity  of  the  memory  storage  of  the  computer  available, 
the  patitioning  method  is  suggested.  However  if  the  members  of  a 
structure  do  not  meet  orthogonally  the  analysis  of  stresses  will  become 
more  complicated. 
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4PPFM0IX 

F«Rr-3  SOURCE  PROGRAM. 

C   C   PROGRAM  1    LFOU   U/5/66 

C      CAL.  CF  FIXED-FND  MOMENTS  OF  A  CURVED-BEAM  UNDER  UNIF.  LOAD 

PI  =  1.  14159265^6 

X=PI/9. 

GK=14.2 

A=X*(GK+1. ) 

B=(GK-1 .  )*SINF(X)»MNF(X) 

C=(GK-1 . )*SINF(X)*COSF(X) 

n=A+CK»X*<COSF(X)+l . )-B*X/ 2.-2. «(GK+1.)*S INFIX )-C 

F=(A+n*X/2.-B+GK*X»5INF(X)-2.»(GK+l. )* < 1 .-COSFI X > ) 

0=1 .78 

1  =  1 

R«200./PI 
10  BN=-0»R**?*(D*< A+C)+B*E> 

TN=fQ*R**2*( E*( A-C)+B*D) 

COA*»2-B**2-C*»2 

BMF=BN/CD 

TMF=TN/CD 

15  FORMAT </6X»lHA.llX.lH8.11X.lHC»llX.lHDtllX«lHE> 
PUNCH  1 5 

16  FCPMAT(5F12.fl) 
PUNCH  16.A.B.C.D.F 

17  F0RMAT(//5X.14HBFNDING-M0MFNT.7X.14HT0RSI0N-M0KENT.13X.J"R/ 1 
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PUNCH  18,BMF,TMF,R 

IF( 1-7)70,30,30 
?0  R=R-7. 

1  =  1  +  1 

GO  TC  10 
JO  STOP 

END 


C   C   PROGRAM  2  LECU   17/4/66 

C      CAl.  OF  STIFFNFSS  MATRIX  INV.(A)  IN  EQ.I9) 

DTMFNMON  A(in,lO)»S(10)»TriO)»H(10tlO) 

INDFXal 

Pt=3.14159?6536 

X=PI/9. 

E  I  =  3??i»000.*863370./ 144000. 

R=200./PI 

READ.N 
5  0=R/(4.*EI ) 

SR=R»0 

rR=R»?R 

GK«14.2 

x?=7.*x 

S?=S!NF(X2) 
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?.]=S  INFIX) 

ci=cosf(x) 

ft  (1  ,1  )=rR«(X?-S2+GK*(6.*X+S2-8.»<n  )l 

A  (1  »?)».e;R#(2.»SS+2.*GK»<2,-2.#Cl-S5)  ) 

A (2.1 )=A( 1,2) 

A(1,3)=SR*(S2-X2+GK*<4.*S1-X2-S2)  ) 

A  1 3  » 1 )  »  A  (  1  ,  3  ) 

A(2.2)=Q»(X2+S2+GK*(X2-S2>  ) 

A(2»3)=Q*<2.*(GK-1.)*SS) 

A(3,?)=A<2,3) 

A(3,3)=0*<X2-S2+GK*(X2+S2) ) 

NN=N-1 

DC  6  1=1 .3 

no  6  j=i  .•» 

6  H(I ,J)=A( I ,J) 
H(l,l  )=1./H( 1,1  ) 
DC  110  M=] ,NN 
<=M+1 
50  DO  60  1  =  1, M 
S( I )=0 
DC  60  J=l ,M 

60  S(  I  )=S(I)+H( I,J)»H( J,K) 

D  =  0 

DC  70  1=1, M 

70  D=D+H(K,r )»S{ I ) 

D=-D+H(K,K) 
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H  (K  •*"  I  =1  ./D 

nc  pr>   1  =  1  ,m 
BO  HII ,<)=-S( I )*H(K*K) 

nc  on  j=i  ,m 

T(J)=0 
DC  90  1  =  1  ,M 
90  T(J)=T( J)+H(K.I )*H( I ,J) 
DC  100  J=1,M 
100  H(K,J)=-T< J)*H(K»K) 
DC  110  1=1, M 
DC  110  J=1,M 
110  H( I ,J)=H( I ,J)-S( I )»H(K,J) 
15  FCRMATI/5X.34HSTIFFNESS-MATRIX  I NV . ( A )  IN  EO.  9. . 19X, 1HR// ) 
PUNCH  15 
160  PUNCH  8,R 

P  rCRMAT(*?X,F1?.6> 
170  PUNCH  9,( (H( I, J) ,J=1,N) ,I=1,N) 
9  FCRMATI3F16.4) 

IFUNDFX-21  25,30,30 
25  INDFX=INDEX+1 
R=R-7. 
GO  TC  5 
30  STOP 
END 

DATA  CF  PROGRAM  2 
3 
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C   C   PROGRAM  3    LFOU   12/4/66 

C      CAL.  OF  STIFFNFSS  MATRIX  IN  EQ.I12) 

DIMENSION  R<3.3).AINV(3.3),C(3.3).BQ(3,3) 
INDFX»1 

P!=3. 1415926536 
X=PI/9. 
R=?00./PI 
IP]  RFAD  1 . INA, IMA, 1MB 

1  FORMATCHfi) 
10?  RFAD  48>< (AINVI I»JI .JM.IMAI .  I  *  1 » I NA ) 
48  FCRMATI3F16.4) 
5  S1=SINF(X) 
Cl=COSF(X) 

an  ,n=-i.o 

B(l,2>=-0. 
B(l»3>»-0. 
R(2.1 )=-R*Sl 
B(?t?)*-Cl 

B(2.?)=+S1 
B(3.1 >=-R»( 1.-C1) 
B(3,?)=-S1 
B(3.3)=-C1 
DO  121  1=1.3 
DO  121  J«l»3 
121  BO( I.J)=B( I.J) 
DC  7  1=1. INA 
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DC  7  J=l . 1MB 

SUM=0. 

DC  6  <=1.IMA 

6  SUMkRQI I »K)*AINV(K.J)+SUM 
C( I  ,J)=SUM 

7  CCNTINUF 

15  FORMAT(/5X,27HSTIFFNESS-MATRIX  IN  EO.  12.  .  19X .  1HR// ) 
PUNCH  1? 
160  PUNCH  8,R 

fl  FCRMATI52X.F12.6) 
170  PUNCH  9-( (C( I .J) »J=1 tIMB)tl«l»INA) 
9  FCRMATI3F16.4) 

IFIINDEX-2)  25.30,30 
25  INDEX=INDEX+1 
R=R-7. 
GC  TC  102 
30  STOP 

F  >in 


DATA  OF  PROGRAM  3 
3        3        3 

20449.9000  -227054.2800  14197.4900 

-227054.2800  3266355.0000  -278175.0500 

14197.4900  -278175.0500  92871.4480 

29003.9020  -286620.0200  17922.0810 

-286620.0300  3669879.6000  -312540.7200 

17922.0ol0  -312540.7200  104344.7600 


C  C       PROGRAM  4    LPOU   12/4/66 

C      CAL.  OF  STIFFNFSS  MATRIX  IN  FO.I16I 

DTMFNStCN  Alio. 10). S(10), TITO). H( 10,10) 
!NDFX=1 

PT=3.1415Q26536 

X=PI/9. 

Et=3?2i.000.»863370./ 144000. 

R=?0O./PI 

RFAD.N 
5  Q=R/(4.»EI ) 

SR=R*0 

CR=R»SR 

GK=14.2 

X2  =  2.*-X 

S?«SINF<X5>) 

S1«SINF(X) 

SS=S1**2 

C1=C05F(X) 

A ( 1 .1 ) =CR*(X2-S2+GK*(6.«X+S2-8.*S1 ) ) 

A (1 .2)=-SR*(2.*SS+2.*GK*(2.-2.*Cl-SS> ) 

A(1.3)=SR*(S2-X2+GK»(4.*S1-X2-S2) ) 

A (2.1 )=A( 1 ,2) 

AI3.1 )=A( 1.3) 

A(2.2)=0»(X2+S2+GK*(X2-S2) ) 

A(2.3)=-Q*(2.*(GK-1. I*SS) 

A(3.?)=A(2,3) 

A(3.3)=0*(X2-S2+GK»(X2+S2) ) 
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NN«N-] 

DC  6  1=1 >* 

PS  6  J  =  l  .•> 
6  H(TtJ)»»UtJI 

Mil ,1 )=1./H( 1 .1 ) 

PC  110  M*l ,NN 

K=M  +  1 
50  Of  60  1=1 ,M 

S(I 1=0 

DC  60  J=1,M 
60  S( I )=S( II+HI I,J)*H(J»K) 

0  =  0 

DC  70  1=1 ,M 
70  D=D+H(K.I )«SI I ) 

n  =  -h+H(if  »(C) 

H(K.if  )  =  1./D 

HC  PO  1=1 »M 
80  H(l  »K)=-SU)»H(K»K) 

DC  00  J=l ,M 

T(J)=0 

DC  90  1=1, M 
90  T( J)=T(J)+H(K,I )*H( I .Jl 

DC  100  J=1,M 
100  H(K.J)=-T( J)*H(K»K) 

DC  110  1  =  1  ,M 

DC  110  J=1,M 
110  H(l ,J)=H( I,J)-S( I )»H(K,J) 
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15  F0RMATI/5X.34HSTIFFNESS-MATRIX  INV.(A)  IN  EO.  9. » 19X, 1HR// ) 
PUNCH  1 5 
16(1  PUNCH  8,9 

8  FWMATI52X.F12.6) 
170  PUNCH  °,< (H( T,J),J=] ,N) .1=1. N) 
O  F0RMATCF16.4) 

IFIINDFX-2)  25.30,^0 
25  INDEX=INDEX+1 
R=R-7. 
GO  TC  5 
30  STOP 
END 

DATA  OF  PROGRAM  4 
3 


C   C   PROGRAM  5    LFOU   12/4/66 

C      CAL.  OF  STIFFNESS  MATRIX  IN  EQ.  (18) 

DIMENSION  B(3,3),AINV(3,3>,C(3»3),BQ(3,3) 

INDEX=1 

PI=3. 1415926536 

X=PI/9. 

R=200./PI 

101  READ  1 ,INA,IMA,IMB 
1  F0RMATI3I8) 

102  READ  48, ( (AINVI I.J) • J=l . IMA ) ♦ I -1 , I NA ) 


- 
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48    FCRMATOF16.4) 

5  S1=SINF(X ) 
C1=CCSF(X) 
B(]  .1 )=-1.0 
B(  1  ,?)=-0. 
B(l  ,3)=-0. 

p(2.?)=-ri 

B(2,3)=-S1 
B(3,1)=-R*(1.-C1) 
B(3.2)=+S1 
B(3.3)--C1 
DC  1?]  1=1,3 
DC  171  Jil.3 
121  RO( I ,J1=B( I ,J) 
DC  7  1=1, INA 
DC  7  J=l ,1MB 
SUM=0. 
DC  6  K=l .IMA 

6  SUM=BO( I,K)*AINV<K,J)+SUM 
C( I ,J)=SUM 

7  CCNTINUE 

15  FCRMAT(/5X,27HSTIFFNESS-MATRIX  IN  EQ.  12. , 19X, 1HR// ) 
PUNCH  15 
160  PUNCH  8,R 

8  FCRMAT(52X,F12.6) 

170  PUNCH  9,1 (C( I»J)»J=1 , 1MB). 1=1, INA) 
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0    F0RMATC»F16.4> 

IFHNOFX-2)     25.30,30 
T5    INDEX=!NDFX+1 

R=R-7. 

GO  TC  102 
30  STOP 

END 

DATA  OF  PROGRAM  5 
3        3        3 

20449. 9f  00 
227056.2600 

14197.4900 

29003.9020 
286620.0300 

17922.0810 


227054.2800 
3266355.0000 

278175.0500 

286620.0200 
3669879.6000 

312540.7200 


14197.4900 

278175.0500 

9287] .4480 

17922.0810 

312540.7200 

104344.7600 


C   C   PROGRAM  6    LFCU   11/26/66 

C      CAL.  OF  STIFFNESS  MATRICES  CF  STRAIGHT  BEAM. IN  E0.(59)»(61) 

DIMFNSICN  FI3i3).G(3>3l 

E I =3220000. #228750./ 144000. 

SL  =  7. 

GJ=322UOOO.*47000./(2,25*144000.) 

F(l,l )=12.*EI/(SL*»3) 

F( ] ,?)=6.*EI/(SL»*2> 

F(l,3)=0. 

F(2»l )=F< 1.2) 

F(2,2)=4.«EI/SL 

F(2,3)=0. 
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FI3.1 1=0. 

F(3.?)=0. 

F(3,3)=GJ/SL 

Gllil )=-12.*EI/ISL**3> 

Gil  .?)=-6.*EI/ISL**2l 

Gl  1.3  1—0. 

G<?,l)«-G(1.2) 

G(?.?)=+2.»EI/SL 

G(2.3>=-0. 

GI3.1 )=-0. 

GI3.?l=-0. 

G(3.3)=-GJ/SL 
1  >  FCRMATI/10X.29HSTIFFNESS-MATRIX  FOR  NEAR  END/ I 

PUNCH  10 
16  FCRMATI3F16.4) 

PUNCH  16.1  IFI  I .J) ,J=1,3) .1  =  1.3) 
21  FCRMAT(//9X,28HSTIFFNESS-MATRIX  FCR  FAR  END/) 

PUNCH  21 

PUNCH  16.1 (Gl I .J) .J=] .3) .1=1.3) 

STOP 

END 


^- 
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C   C   PROGRAM  7    UFCU   11/26/66 

C      CAL.  OF  STIFFNESS  MATRICFS  OF  STRAIGHT  BEAM.  IN  EQ. (66)1(68) 

DIMFNSION  F(  3.3)  .GO. 3) 

F I =3220000. *228750./ 144000. 

SL  =  7. 

GJ=322COOO.*47000./(2.25*144000.) 

Fll.l  )=12.*E!/(SL*»3> 

F(]  .2)=-6.*EI/(SL*»2> 

F(l  .31=0. 

F(?»1 )*F( 1.2) 
F(?,?)=4.*EI/SL 
F(2.3)=0. 
F(3.1)«0. 
F(3,2)=0. 
F(3.3)=GJ/SL 
G( 1.1 )=-12.»EI/(SL»*3> 
G(l  .2)=6.*EI/(SL**2) 
G(  l,3)=-0. 
GI2.1 )=-G(l .2) 
G(2.?)=+2.*EI/SL 
G(2.3)=-0. 
GI3.1 )=-0. 
Gt3,?)=-0. 
G(3.3)=-GJ/SL 
10  F0RMAT(/10X,29HSTIFFNESS-MATRIX  FOR  NEAR  END/1 
PUNCH  10 
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16  F0RMATI3F16.4) 

PUNCH  16iMFUtJ)»J«l»3).f«l»3) 

21  FCRMAT(//9X,28HSTIFFNESS-MATRIX  FOR  FAR  END/) 
PUNCH  21 

PUNCH  16,< (G( I,J) ,J=1,3) ,1=1,3) 
STOP 
END 


C   C   PROGRAM  8   MATRIX  MULTIPLICATION    LEOU   11/25/66 
DIMENSION  A(12.?M.3l?4,12l.fllM2l 
100  READ  1 ,INA, IMA, I  MR 

1  FORMAT(3I8) 

READ  4,( (A( I ,J) ,J= 1,1  MA), 1  =  1, INA) 
READ  2,{ (R( I ,J) ,J=1, 1MB) ,1=1, IMA) 
h    F0RMATI6F11.8) 

2  F0RMATI6F12.2) 
DO  7  1=1 , INA 
DO  7  J=l ,IMB 
SUM=n. 

10  6  tf  =  l  ,IMA 
ft  SUM=« ( I ,K)*P(K,J)+SUM 

C(I ,J)=SUM 
7  CONTINUE 

PUNCH  3,(  (C(  I,J),J  =  1 , I  MP), 1=1, INA) 

3  FCRMATI6F12.2) 
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STOP 
END 
C      DATA  OF  PROGRAM  8 
1?       24       12 


C   r   PROGRAM  9   MATRIX  INVERSION  RY  PARTITIONING    LFOU   11/2/66 

DIMENSION  A(20.20),B( 19) ,C<19) 
4  FORMATI6F12.?) 

RFAD.N 

NN=N-1 
10  RFAD  4,( ( A< I ,J) ,J=1,N> ,I=1,N) 

A( 1,1 )=1./A< 1.1) 

DO  110  M=1,NN 

K«M+1 
50  DO  60  1=1, M 

B< I  1=0 

no  60  j=i  ,m 

60  B( I )=B< I )+A( I,J)*A(J»K) 

0  =  0 

DO  70  1  =  ]  ,M 
70  D  =  D  +  A(K, I )»B(  I  ) 

D=-D+A (K,K> 

MK»K)»1./D 

DO  80  1  =  1, M 
80  All  ,K)=-B( I )*A(<,K) 
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dc  oo  j=i (m 

C( J)=0 
DC  90  I=1,M 
00  C(J)=C< J)+AIK,I )*A( I ,J> 
DC  100  J=l |M 
100  A(K»J)=-C( J)*A(K»K) 
DC  IK'  1=1. M 
DC  110  J»1»M 
110  A< I  .J)=A(  I .J)-R( I )*A(K.J) 

DC  170  1=1, N 
170  PUNCH  ^ , ( A ( I , J ) ,J=1 ,N) 
5  FCRMATUE16.8) 
END 
C      DATA  CF  PROGRAM  9 
1? 


C   C   PRCGRAM  10   MATRIX  TIMES  VECTOR    LECU   11/29/66 
DIMENSION  A  I  48,1?) ,VEC(12) .PRCDU8) 
10O  rfao  1 .IMA. IMA 

1  FCRMATI2I8) 

2  FCRMATI6F12.2) 
?  FCRMATIF20.15) 

RFAD    2.( (A( I.J) .J=l . IMA). 1=1. INA) 
READ    1.IVECI I ) .1=1. IMA) 
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DO  115  1=] »!NA 
SUM=O.U 

DO  lift  <  =  1 » I  MA 
116  SUM  =  A(  I,<)*VEC(ICI  +  SUM 

ppodi i i=sum 

115  CCNTINUF 

5  FrRMAT(?nx,3?HPR0DUCT  OF  A  MATRIX  AND  A  VECTOR//) 

PUNCH  5 

PUNCH  6. (PRODI  I ) ,1  =  1 ,INA| 
ft  FORM\T(30X.F10.5) 

STOP 

END 

DATA  OF  PROGRAM  10 

48        12 


C   C   PROGRAM  1]    LEOU   11/25/66 

C      CAL.  OF  CONSTANTS  G.H.S  IN  E0.(46).(47)  AND  (48) 

PI=3. 1415926536 

IMDFX=1 

X=PI/9. 

P*l. 

GK=14.2 

R=200./PI 
5  Z=PI/36. 

DZ»PI/36. 
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SX=SINF (X) 

X?=2.*X 

SX?»S INFIX?) 
CX-CCSFIX) 

r x?=ro^P( x?i 

y=gk+i . 

YN=Gkr-l. 
10  SZ=SINF(Z) 
Z?=2.*Z 
CZ=CCSF(Z) 
SZ2=SINF<Z2) 
CZ2=CCSF(Z2> 

xz=x-z 

SXZ=SINF(X-Z> 

E I = 3220000. *863370./ 144000. 

T=(P»R«»2)/(4.*FI ) 

CP=P»T 

?o  r,*  =  -YN»rz*f;x2-sz?) 

GR=SZ*YN* ( CX2-CZ  2 ) 

C,r  =  -?.»XZ*(  Y*CZ+2.*GK  ) 

G0=4.*GK*SXZ+4.*GK»{ SX-SZ ) 

G» (GA+GB+GC+GD)*CR 

HA=2.*XZ»Y*SZ 

HB=4.*GK#(CX-CZ) 

HC=-CZ*(CX2-CZ2)»YN 

HD=-SZ*<SX2-SZ2)*YN 

H=T»(HA+HB+HC+HDI 
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SA=2.*XZ*Y»CZ 
SB«-4.«GK«(SX-SZ) 
SC=CZ*(SX2-SZ2)*YN 
SD=-<;Z»<CX2-CZ2)*YN 
S«T*(SA+Sfl+SC+SO) 
25  FORMATI  nX.1HR.Fl?. 8) 
PUNCH  25, R 

30  'rCRM«T(//4X,lHG,15X.lHH.15X.lHS,19X,lHZ) 

PUNCH  30 
35  FCRMATI50X.F14.8) 

PUNCH  35. Z 
AC  F~RMAT(3F16.12> 

PUNCH  40.G.H.S 

IFIZ-PI/12.)  50,54,54 
50  Z=Z+DZ 

GO  TO  10 

54  CONTINUE 
IHINDFX-2)     55,60,60 

55  !NDFX=INOFX+l 
R=R-7. 

GO    TC    5 
60    STCP 
FND 
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C   C   PROGRAM  1?    LFCU   11/29/66 

C      CAL.  OF  FIXED-FND  STRESSES  IN  FQ.I45) 

DIMENSION  A(12.12).VECI12> .PRODI  12) 

PI  =  3.  1415926536 

DZ»PI/36. 

N  =  l 

r=?oo./pi 

IOC  RFAD  l.INA.IMA 
]  F0RMATI2I8) 
?  F0RMATMF16.4) 

3  F0PMATI3F16.1?  I 
16  Z=PI/36. 

INDEX=1 
B  RFAD  ?.( < A! I,J),J=1.IMA).I=1.INA) 
9  RFAD  3.IVFCI  I ) .1  =  1 » I  MA ) 
DO  115  1=1. INA 
SUM=0.0 

DO  116  K=1.IMA 
116  SUM=A( I .K)»VFC(K)+SUM 
PRODI  I  l=SUM 
'  ll J  rCNTINUF 

4  FORMAT! /I 3X.1HR.F15.8) 
PUNCH  4,9 

5  FORMAT! 17X.18HFIXFD-END  STRESSES) 
PUNCH  5 

4  0  F0PMATI/6X.1HF.13X.1HM.13X.1HT.15X.1HZ) 
PUNCH  40 
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6    F0RMAT(3FI4.8»F1A.8) 

PUNCH    6. (PRODI  I )tI  =  l.INA).Z 

IF( INDEX-3l20t30»30 
?0    IND,iX"INOEX+l 

Z=Z+DZ 

GC    TO    9 
30    IF(N-2)15t50,50 
15    N  =  N  +  1 

P=R-7. 

Gr.  Tr:   16 

SO    STOP 
FND 

DATA    CF    PROGRAM    12 

3  3 

-20449. 90U0  +227054.2800  -14197.4900 

+227054.2800  -326635  5.0000  +278175.0500 

-14197.4900  +278175.0500  -92871.4480 

-.000140844240  -.000009378848  -.000006671969 

-.000066033581  -.000003845325  -.000001600150 

-.000017110796  -.000000864254  -.000000061282 

-29003.9020  +286620.0200  -17922.0810 

+286620.0300  -3669879.6000  +312540.7200 

-17922.^810  +312540.7200  -104344.7600 

-.000099305638  -.000007429725  -.000005285392 

-.000046558573  -.000003046185  -.000001267605 

-.000012064381  -.000000684644  -.000000048  546 
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C   C   PROGRAM  n    LEOU   12/7/66 

C      CAL.  OF  FIXED-END  STRESSES  IN  EO.(bO) 

DIMENSION  B(3.3>»AINV<3»3>»C(3.3).t)QI3«3).D13.3>tFI3>.G<3>»H(3) 

DIMENSION  VEC(3) .PRODI 3) 

INDEX=1 

P  =  l. 

PI=3. 1415926536 

X  =  PI/<3. 

R=20o./pi 

101  rfad  1 < tna, ima. 1mb 

1  F0RMATI3I8) 
?6  2=PI/36. 
DZ=PI/36. 

102  RFAD  48, ( (AINVI I »J) .J=1.IMA> .1=1. INA) 
48  FCRMATI3F16.4) 

5  S1=SINF(X) 
C1=C0SF(X) 
XZ=X-Z 

SXZ=SINF(XZ> 
CXZ=rCSF(XZ) 
B(  1.1  )=+1.0 
B(l.?)=+0. 
B(1.3l=+0. 
BI2.1 )=+R*S] 
B(2,2)=+C1 
B(2.3)=-S1 
B(3,1)=+R»(1.-C1I 
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R(3»?l=+S] 
B(3»3)»+C1 

00  171  1=1 ,3 
DC  1?1  J=1.3 
1?1  RO( I.J)=f( I .J) 
DC  7  1  =  1 ,  INA 
DC  7  J  =  l tIMR 
SUM=0. 
DC  6  K  =  l .IMA 

6  SUM=BQ( I ,K)*AINV(K.J)+SUM 
C( I ,J)=SUM 

7  CONTINUE 

DC  130  1=1,3 

DC  130  J=l,3 
130  D( I  ,J)=C(  I ,J) 
10  READ  15, ( VECI I ) ,1=1, IMA) 
15  FCRMAT(*F16.1?> 

DC  115  1=1 , INA 

SUM=D.O 

DC  116  K=l .IMA 
116  SUM=D( I ,K)*VEC(K)+SUM 

PPCDI I )=SUM 
115  CONTINUE 

DO  120  1=1,3 
120  F( I )=PRCD( I ) 

G(1I=P 

G(2)=P»R*SXZ 
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G(3!=P*R«(1.-CXZ> 

no  115   1=1.3 

135    HI  I >=F( I )+G(  I ) 

16  F0PMATI/13X.1HR.F15.8) 
PUNCH    16.  P. 

17  F0RMATU7X.18HFIXFD-FND    STRFSSFS) 
PUNCH    17 

41    F0RMATI/6X.1HF.13X.1HM.13X.1HT.15X.1HZ) 

PUNCH    41 
20    F0RMATI3F14.8.F14.8) 

PUNCH    Z0.(H(I)»1»1»INA»»Z 

IFIZ-PI/12.)     30.40.40 
30    Z=Z+DZ 

GO    TO    5 
40    CONTINUF 

IFUNDFX-2)     45.60.60 
46    INDFX=IN0FX+1 

RsR-7. 

GO  TO  26 
60  STOP 

FND 

DATA  OF  PROGRAM  13 

3        3  3 

20449.9000  -227054.2800  14197.4900 

-227054.2800  3266355.0000  -278175.0500 

14197.4900  -278175.0500  92871.4480 

-.000140844240  -.0001/09378848  -.O0OOU6671969 

-.000066033581  -.000003845325  -.000001600150 

-.000017110706  -.000000864254  -.000000061 2R2 

29003.9020  -286620.0200  17922.0810 

-286620. C300  3669879.6000  -312540.7200 

17922.0810  -312540.7200  104344.7600 

-.000099305638  -.000007429725  -.0QO0U52R5392 
-.000046558573  -.000003046185  -.000001567605 
-.000012064381   -.000000684644   -.000000048546 


Table  1 .  Final  member  end  stresses  of  the  structure  in  Fig.  8  due  to 
uniform  load. 

fixed 
j  P.  I  +  end   -  final  member  end  stresses 
'   stresses 

F^  47.65  +  19.778  -  +  67.428  kips 

k\  -767.50  -  75.369  -  -  842.829  ft-kips 

t}  56.73  -  O.347  =  +  56.383  ft-kips 

F1 

2  -47.65  +  19.778  =  -  27.872  kips 

Mj  -296.85  +  75.369  -  -  221.481  ft-kips 

T^  26.26  -  O.347  -  +  25.913  ft-kips 

?2  3.75  +  4.340  .  +   8.090  kips 

M^  20.64  +  5.060  .  +  25.7OO  ft-kips 

T£  6.10  +  0.000  -  +   6.100  ft-kips 

F2  -  .00  +  19.778  -  +  19.778  kips 

Mg  290.75  -  75.369  -  +  215. 381  ft-kips 

?22  -51.27  -  0.347  -  -  51.617  ft-kips 

Fj  .00  +  19.778  -  +  19.778  kips 

Mj  -290.75  +  75.370  -  -  215.381  ft-kips 

T^  -51.27  -  O.347  .  -  51.617  ft-kips 

Fj  3.75  +  4.340  -  +   8.090  kips 

K 

K,  20.64  +  5.060  =  +  25.700  ft-kips 

T^  -6.10  +  0.000  =  -   6.100  ft-kips 

F^  -47.65  +  19-778  -  -  27.872  kips 

M^  296.85  -  75.369  ■  +  221.481  ft-kips 

t|  26.26  -  O.347  =  +  25.913  ft-kips 

F^  47.65  +  19.78O  -  +  67.428  kips 

M^  767.50  +  75.369  -  +  842.829  ft-kips 

t'  56.73  -  O.347  -  +  56.383  ft-kips 


DC 


Table  2.  Final  member  end  stresses  of  the  structure  in  Pig.  8  due  to 
uniform  load. 


fixed 

end 

Btresses 


JP_l  +  end   ■  final  member  end  stresses 


F5  35.79  +  17.605  =  +  53.403  kips 

M^  -543.80  -  59.706  .  -  603.536  ft-kips 

T^  34.56  -  0.275  -  +  34.265  ft-kips 

Fg  -35.79  +  17.603  -  -  18.197  kips 

Mg  -170.81  +  59.706  -  -  111.144  ft-kips 

Tg  31.21  -  0.275  '   +  30.925  ft-kips 

Fg  -3-75  +  4.340  =  +   O.57O  kips 

Mg  5.59  -  5.060  .  +   O.55O  ft-kips 

Tg1  -6.10  +  0.000  =  -   6.100  ft-kips 

Fg"  -  .00  +  17.603  -  +  17.603  kips 

H^  176.92  -  59.706  -  +  117.244  ft-kips 

Tg  -31.19-  0.275  -  -  31.475  ft-kips 

pI  .00  +  17.603  >  +  17.603  kips 

M^  -176.92  +  59.706  =  -  117.244  ft-kips 

T7,  -31.19  -  0.275  =  -  31.475  ft-kips 

Fy  -3.75  +  4.340  -  +   0.570  kips 

m5  5.59  -  5.O6O  =  +   O.55O  ft-kips 

T:j  6.10  +  0.000  =  +   6.100  ft-kips 

Fy  -35.79  +  17.603  ~  -  18.197  kips 

M^  170. 81  -  59.706  -  +  111.144  ft-kips 

T?,  31.21  -  0.275  -  +  30.925  ft-kips 

Fg  35.79  +  17.603  .  +  53.403  kips 

Mg  543.80  +  59.706  .  +  6O3.536  ft-kips 

Tg*  34.56  -  0.275  -  +  34.285  ft-kips 
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TAPl.F  *,  MATRIX  0 

--UNIT  1  CAD  AT  POINTS  9,10,11, 

j 

UNIT 

.CAD  AT  POINT 

Q 

10 

1  I 

7 

F? 

-. 19453180 

-.50000000 

-.84544850 

-1.00000000 

VU 

-1.07909800 

-2.8781 1500 

-3.20O24310 

0.00000000 

T2 

.00822890 

.01644997 

.01027509 

0.00000000 

F3 

0.OCO00O00 

O.OOOOOOOO 

O.OOwOUOOO 

0.00000000 

H3 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

T3 

o.oooooooo 

0.00000000 

0.00000000 

0.00000000 

F6 

o.pooonooo 

0,00000000 

0.00000000 

0.00000000 

M6 

0,00000000 

O.uCOOOOOO 

o.oooooooo 

0.00000000 

?6 

O.OOOOOuOO 

O.OOOOOOOO 

O.OOOOwOOO 

0.00000000 

F7 

o.ooooouou 

0.00000000 

0.00000000 

0.00000000 

M7 

0.00000000 

0.00000000 

o.oooooooo 

0.00000000 

t? 

0 .oonnnoon 

0.00000000 

O.oooooooo 

0,00000000 

table  4.  t>,ATRix  q — unit  load  at  points   12,13,14,3 


UNIT 

LOAD  AT  POINT 

1? 

13 

14 

3 

p2 

-.84546821 

-.50000000 

-.15465155 

0.00000000 

M2 

3.20061630 

2.87812000 

1.07916620 

0.00000000 

T2 

.01030862 

.01644585 

.00820781 

0.00000000 

F3 

-.15453180 

-.50000000 

-.84544850 

-1.00000000 

M5 

-1.07909800 

-2.87812000 

-3.20024310 

0.00000000 

T3 

.C0822890 

.01644997 

.01C27509 

0.00000000 

F6 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

M6 

0.00000000 

O.OuOOOOOO 

O.OOOOOOOO 

0.00000000 

T6 

0.00000000 

O.OOOOOOOO 

0.00000000 

0.00000000 

F7 

0,00000000 

0.00000000 

0.00000000 

0.00000000 

M? 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

T7 

0.00000000 

0.00000000 

0.00000000 

0.00000000 
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TSPLE  5.  MATRIX  C 

— UNIT  LOAD  AT 

POINTS  15.16.17 

.18 

UNIT 

LCAD  AT  POINT 

15 

16 

17 

18 

P2 

0.00000000 

O.OUOOuOOO 

O.OOl/UOuOO 

-.50000000 

M2 

0.00000000 

o.oooooooo 

0.00000000 

-.87500000 

T2 

o.rnoooooo 

o.oconnooo 

0.00000000 

o.oooooooo 

p3 

-.84546821 

-.50000000 

-.15455155 

0.00000000 

K, 

3.20061630 

2.87812000 

1 .07916620 

0.00000000 

'h 

.01030862 

.01644585 

.00820781 

0.00000000 

F6 

0.00000000 

0.00000000 

0.00000000 

-.50000000 

M6 

0.00000000 

0.00000000 

0.00000000 

.87500000 

T6 

o.ooonAoio 

o.oooooooo 

0.00000000 

o.oooooooo 

P7 

o.of  o^nooo 

0.00000000 

0,00000000 

o.oooooooo 

M? 

n. oooooooo 

o.oooooooo 

o.oooooooo 

0.00000000 

T7 

0.00000000 

0.00000000 

O.OOUOOOOO 

0.00000000 

TaBLE  6.  MATRIX  0 — UNIT  LOAD  AT  POINTS  19.20.21 .22 


3 


M7 
In 


19 


o«ooooococ 

0,00000000 

o.ooocccco 

-.50000000 
-.87500000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 

-.50000000 

.87500000 

0.00000000 


UNIT  LOAD  AT  POINT 


20 


0.00000000 
0.00000000 
o.ocooooco 
o.oooouooo 

0.00000000 

0.00000000 

-.15453170 
-.96044500 
.00732360 
0.00000000 
0.00000000 
0.00000000 


21 


0.001C00O0 
O.OOCOOOOO 
O.OOOOOOOO 
O.OOwuUOOO 
0.00000000 
0.00000000 

-.50000000 

-2.56165000 

.01464100 

0.00000000 

O.OOUUOOOO 

0.00000000 


22 


0.10000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
-.84544860 
-2.84835840 
.00914525 
0.00000000 
0.00000000 
0.00000000 
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T*FLF  7.  MATRIX  0 — UNIT  | CAD  AT  POINTS  6*23.24.25 


UNIT 

LOAD  AT  POINT 

6 

23 

24 

25 

F2 

o.opr.opnoo 

0.00000000 

o.oococooo 

o.oooooooo 

M2 

n.noonroon 

O.OPOOPOOO 

0.00000000 

o.ooonoooo 

?2 

o.ppopoooo 

0.00000000 

0.00000000 

0.00000000 

P3 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

M* 

0.00000000 

O.CCOOOOOO 

0.00000000 

0.00000000 

T3 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

P6 

-1 .ooPPOPPO 

-.84546832 

-.50000000 

-.15455148 

H6 

P.P0POO0P0 

2.84868920 

2.56165000 

.96050480 

T6 

O.OOPPOOPO 

.00917510 

.01463761 

.00730520 

P7 

0.00000000 

-.15453170 

-.50000000 

-.84544860 

M7 

0.00000000 

-.96044500 

-2.56165000 

-2.84835840 

T7 

0.  10000000 

•0073Z36C 

.01464100 

.00914525 

TSF 

LF  8.  MATRIX  0 

--UNIT  LOAD  AT 

POINTS  7,26,27, 

?fi 

UNIT 

_OAD  AT  POINT 

7 

26 

27 

28 

*2 

0.00000000 

1   0.00000000 

0.00000000 

0.00000000 

M2 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

T2 

O.OOOOOOOO 

0.00000000 

0.00000000 

0.00000000 

P3 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

M* 

o.ppoonooo 

0.00000000 

0.00000000 

0.00000000 

h 

0,00000000 

o.oooooooo 

o.ooocoooo 

0.00000000 

p6 

P.nORPOOOO 

0.00000000 

0.00000000 

0.00000000 

M6 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

T6 

0.00000000 

o.oooooooo 

0.00000000 

0.00000000 

F7 

-1.00000000 

-.84546832 

-.50000000 

-.15455148 

M7 

o.oooooooo 

2.84868920 

2.56165000 

.96050480 

T7 

o.oooooooo 

.00917510 

.01463761 

.00730520 
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T^PtF    9.    FINftL    MFMRfR    fno    STRFSSES 

UNIT    LOAD    AT    POINTS    9«10>1 1 •?• 12i 


13 


UNIT  LOAD 

AT  POINT 

0 

10 

11 

2 

12 

13 

pl 

.97255 

.89568 

.78215 

.64993 

.51628 

.38767 

Ml 

-4.74)57 

-7.83342 

-9.38622 

-9.69314 

-9.13661 

-7.90639 

f\ 

,05070 

.19010 

.36607 

.49988 

.54712 

.52594 

»1 

.02745 

. 10432 

.21785 

.35008 

-.51628 

-.38767 

1  2 

-.21906 

-1.01026 

-2.53080 

-4.87174 

-2.46867 

-.83152 

«1 

-.02643 

-.03488 

.06966 

.35026 

.62862 

.72155 

»s 

-.01223 

-.04269 

-.07975 

-.10904 

-.11888 

-.11233 

m| 

.01797 

.04569 

.05203 

•0276e 

-.01683 

-.08219 

if 

-.00299 

.00009 

.02681 

.09277 

.15780 

.17696 

»i 

-.01522 

-.06162 

-.13810 

-.24104 

.63516 

.50000 

M| 

.22205 

1.01018 

2.50399 

4.77897 

2.310R7 

.65457 

t| 

.00846 

-.01081 

-.12168 

-.37793 

-.61178 

-.63936 

F| 

.01522 

.06162 

.13810 

.24104 

•364R4 

.50000 

M? 

.12558 

.38881 

.61230 

.62823 

.27344 

-.65457 

t| 

-.02547 

-.09875 

-.21187 

-.35396 

-.52064 

-.63936 

P5 

-.00565 

-.02107 

-.04339 

-.06942 

-.09436 

-.11233 

M3 

-.02362 

-.07759 

-.13472 

-.16840 

-.14626 

-.08218 

T5 

-.01124 

-.04070 

-.08081 

-.12354 

-.16056 

-.17696 

Fi 

-.00957 

-.04055 

-.09471 

-.17161 

-.27048 

-.38767 

M3 

-.11434 

-.34812 

-.53149 

-.50469 

-.11288 

.83152 

t2 

,04°09 

. 17634 

.34660 

.52236 

.66691 

.72155 

pi 

.00957 

.04055 

.09471 

.17161 

.27048 

.38767 

M| 

.33252 

1.27032 

2.68008 

4.38952 

6.22349 

7.90638 

Tl 

.02970 

.10903 

.21969 

.34062 

.45036 

.52594 
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r«MTlM(j/iTi?M  "F  TABl.F  9 

UNIT  LOAr 

AT  POINT 

o 

10 

1 1 

2 

12 

13 

p6 

.01?34 

.04271 

.07911 

.10757 

.11 772 

.11233 

VL% 

-.29205 

-1.00995 

-1 .87921 

-2.59851 

-2.93084 

-2.87971 

T6 

-.01240 

-.03935 

-.06534 

-.07698 

-.06906 

-.04656 

F6 

-.01234 

-.04271 

-.07911 

-.10757 

-.11772 

-.11233 

b 

.03098 

.10797 

.21039 

.33076 

.44903 

.51323 

m6 

v6 

.06936 

.23647 

.43380 

.59349 

.66503 

.64482 

F£ 

.01223 

.04269 

.07975 

.10904 

.11888 

.11233 

4 

-.10359 

-.34453 

-.61028 

-.79093 

-.81530 

-.70412 

?l 

.00299 

-.00009 

-.02681 

-.09277 

-.15780 

-.17696 

A 

.00011 

.00001 

-.00064 

-.00146 

-.00115 

-0.00000 

4 

-,03'»96 

-.10788 

-.18358 

-.23799 

-.29123 

-.33627 

n 

.03423 

.10806 

.17648 

.19744 

.15027 

.05930 

y 

-.00011 

-.0000] 

.00064 

.00146 

.00115 

0.00000 

"? 

,04]44 

.13808 

.24527 

.31952 

.34737 

.33628 

T? 

-.02093 

-.06470 

-.10087 

-.09913 

-.03767 

.05930 

P? 

.00565 

.02107 

.04339 

.06942 

.09436 

.11233 

m5 

-.01594 

-.06994 

-.16901 

-.31757 

-.51424 

-.70412 

T^ 

.01124 

.04070 

.08081 

.12354 

.16056 

.17696 

F^ 

-.00554 

-.02106 

-.04403 

-.07089 

-.09551 

-.11233 

M^ 

-.0^268 

-.17877 

-.32608 

-.44306 

-.50793 

-.51323 

m8 

x7 

.03687 

.13463 

.26988 

.41670 

.55191 

.64482 

P§ 

.00554 

.02106 

.04403 

.07089 

.09551 

.11233 

M§ 

.16945 

.62220 

1.25200 

1.93260 

2.51697 

2.87970 

J$. 

.00730 

.00660 

.00838 

.00219 

-.01854 

-.04656 
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TiP|_r    lr.    FF*'AL    "^"PrR    pvn    c,TP!~SSF'"- 

I.WIT    LOAD    AT    POINTS    14« 3 .1 5 t 16 » 17»I 8 


UNIT  LCAT 

AT  POINT 

14 

3 

15 

16 

17 

18 

Fi 

.2TP48 

.17161 

.09470 

.04055 

.00957 

.44215 

mJ 

-6.223*5 

-4.38952 

-2.68006 

-1.27032 

-.33258 

-6.03375 

Tl 

.45036 

.34062 

.21969 

.10903 

.02971 

.39245 

F| 

-.27048 

-.17161 

-.09470 

-.04055 

-.00957 

-.44215 

H| 

.11281 

.50469 

.531 52 

.34812 

.11434 

-3.82314 

Tl 

.66689 

.52236 

.34660 

.17634 

.04910 

-.00266 

Pi 

-."9436 

-.06942 

-.04339 

-.02108 

-.00565 

.56592 

M^ 

-.14628 

-.16840 

-.1 3473 

-.07759 

-.02362 

1 .19219 

T4 

.160^6 

.12354 

.08081 

.04070 

.01124 

.00405 

pi 

.36484 

.24104 

.13810 

.06162 

.01522 

-.12377 

m| 

-.27337 

-.62822 

-.61232 

-.38881 

-.12559 

2.94409 

■I 

-.52061 

-.35396 

-.21187 

-.09875 

-.02548 

-.31453 

'! 

.63516 

-.24104 

-.13810 

-.06162 

-.01522 

.12377 

m| 

-2.31095 

-4.77896 

-2.50391 

-1 .01017 

-.22210 

-.17917 

T2 
X3 

-.61173 

-.37793 

-.12166 

-.01081 

.00846 

-.23618 

P^ 

-.11888 

-.10903 

-.07975 

-.04269 

-.01224 

.01044 

M5 

-.01685 

.02768 

.05206 

.04569 

.01795 

.04690 

T5 
x3 

-.15779 

-.09277 

-.02680 

-.00009 

.00298 

-.02532 

p3 
3 

-.51628 

,35007 

.21785 

.10432 

.02745 

-.13421 

m| 

2.46874 

4.87174 

2.53071 

1 .01026 

.21912 

.20448 

*3 

.62859 

.35026 

.06961 

-.03488 

-.02641 

.18929 

>? 

•  M628 

.64993 

.78215 

.89568 

.97255 

.13421 

m| 

9.13649 

9.69313 

9.38639 

7.83341 

4.74136 

2.79480 

*I 

.54711 

.49988 

.36606 

.19010 

.05073 

.26745 
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r 

"MTrwu&Tir 

H    OF  TAR|_C 

10 

U*IT  LCAO 

AT  POINT 

14 

3 

15 

16 

17 

18 

p6 

.09551 

.07089 

.04403 

.02106 

.00554 

.30251 

m| 

-7.51695 

-1.93260 

-1.25199 

-.62220 

-.16949 

-4.18353 

T6 

-.01853 

.00219 

.00838 

.00660 

.00230 

.27819 

p6 

-.09551 

-.07089 

-.04403 

-.02106 

-.00554 

-.30251 

146 

b 

.50790 

.44306 

.32609 

.17877 

.05269 

-1.83615 

T6 

•*51°0 

.41670 

.26988 

.13463 

.03688 

.13571 

^ 

.09436 

.06942 

.04339 

.02108 

.00565 

.43408 

4 

-.51424 

-.'1757 

-.16907 

-.06994 

-.01594 

-.73074 

< 

-.16056 

-.17354 

-.0808) 

-.04070 

-.01124 

-.00405 

*l 

.00115 

.00146 

.00064 

-.00001 

-.00011 

-.13157 

4 

-.34734 

-.31952 

-.24528 

-.13808 

-.04145 

2.71521 

n 

-.03767 

-.09913 

-.10087 

-.06470 

-.02094 

-.27997 

FT 

-.00115 

-.00146 

-.00064 

.00001 

.00011 

.13157 

mT, 

.29121 

.23800 

.18361 

.10788 

.03395 

-.09752 

tT, 

.l"iC26 

.19744 

.17649 

.10806 

.03423 

-.21598 

* 

.1  1888 

. 1C903 

.07975 

.04269 

.01224 

-.01044 

U5 

-.81 *?7 

-.79093 

-.61030 

-.34453 

-.10360 

.02618 

■? 

.15779 

.09277 

.02680 

.00009 

-.00298 

.02532 

pB 
7 

-.1 1772 

-.10757 

-.07911 

-.04771 

-.01235 

-.1211? 

M§j 

-,440pn 

-."077 

-.21041 

-.10796 

-.03097 

.07770 

T8, 

.66*01 

.59349 

.43380 

.73647 

.06937 

.18981 

*8 

.11 777 

.10757 

.0791.1 

.04271 

.01235 

.17113 

M8 
u8 

2.93087 

2.59851 

1 .87971 

1.00995 

.2921 1 

2.34446 

T§ 

-.06906 

-.07698 

-.06535  ■ 

-.03935 

-.01240 

.21085 
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UN 


IT    LOAD    AT 


i5   rron   ST"17 
POINTS    1" 


20,21  »?2.6>23 


UNIT  LCAD 

AT  POINT 

19 

20 

21 

22 

6 

23 

»1 

.14659 

.02236 

.07521 

.13555 

.18094 

.19691 

Ml 

-3.17143 

-.25630 

-.86872 

-1.58509 

-2.16366 

-2.44253 

Tl 

.27938 

.03364 

.11506 

.21176 

.29050 

.32618 

»i 

-.14659 

-.02236 

-.07521 

-.13555 

-.18094 

-.19691 

M^ 

-.11600 

-.23453 

-.78202 

-1.38944 

-1.80724 

-1.8806) 

■i 

.25937 

-.02980 

-.09977 

-.17726 

-.22766 

-.22710 

p| 

-. 31 K27 

.03057 

.10089 

.17681 

.22553 

.22558 

m| 

-."C39 

.07292 

.24461 

.43692 

.56892 

.58197 

I* 

.04646 

.00640 

.01391 

.00520 

-.03504 

-.07639 

F2 

.16186 

-.00821 

-.02568 

-.04126 

-.04459 

-.02867 

M| 

.06955 

.22813 

.76811 

1.38424 

1.84228 

1.95700 

■i 

-.25598 

-.04311 

-.14484 

-.25966 

-.34126 

-.35487 

*§ 

-.16186 

.00821 

.02568 

.04126 

.04459 

.02867 

h| 

-3.67716 

-.05037 

-.21227 

-.49120 

-.87693 

-1.33601 

■  T2 
x3 

-.40466 

-.00599 

-.02803 

-.07103 

-.13821 

-.22577 

pi 

.55058 

.00039 

.00679 

.02638 

.06459 

.12302 

M5 

1 .1544] 

.00712 

.03791 

.10392 

.21191 

.35575 

T3 

-,05?68 

.00340 

.01372 

.03040 

.05176 

.07367 

Pi 

-.'8872 

-.00860 

-.03247 

-.06764 

-.10918 

-.15169 

H^ 

2.85584 

.04697 

.19855 

.46080 

.82517 

1.26234 

T2 

.12526 

-.00113 

-.00988 

-.03289 

-.07370 

-.12997 

Pl 

.38872 

.00860 

.03247 

.06764 

.10918 

.15169 

m| 

5.82304 

.14267 

.51705 

1.02848 

1.57670 

2.07220 

*2 

.39794 

.01800 

.06604 

.13299 

.20621 

.27277 

95 


r 

"NTtMiiATrr 

N  OF  TA^LF 

1  1 

UNIT  LOAD 

AT  POINT 

l<5 

20 

21 

22 

6 

23 

A 

.10553 

.95837 

.85010 

.70468 

.55410 

.42294 

Mf 

-2.186  tl 

-3.90090 

-5.90906 

-6.46161 

-6.13498 

-5.46658 

T6 

.15678 

.06375 

.23299 

.44213 

.59874 

.64469 

^ 

-. 10553 

.04163 

.14990 

.2953? 

.44589 

-.42294 

M6 

.06305 

-.21382 

-.99286 

-2.48978 

-4.7684* 

-2.83897 

m6 

l6 

.23986 

-.10177 

-.29878 

-.42985 

-.35779 

-.18137 

F£ 

.01 527 

-.03057 

-.10089 

-.17681 

-.22553 

-.22558 

K 

-.10352 

.14108 

.46162 

.80073 

1.00982 

.99710 

Ti 

-.0*646 

-.00640 

-.01391 

-.00520 

.03504 

.07639 

4 

.09026 

-.01106 

-.04901 

-.11851 

-.22036 

.64952 

M7 

-.01659 

.22022 

1.00677 

2.49498 

4.73341 

2.76258 

T? 

-.13634 

-.03930 

-.16284 

-.37088 

-.65204 

-.81573 

1 

-.09026 

.01106 

.04901 

.11851 

.22036 

.35148 

m7 
7 

-1 .78021 

-.00604 

-.05194 

-.17462 

-.40045 

-.77154 

«? 

-.17463 

-.00059 

-.02384 

-.09984 

-.25319 

-.49558 

F? 

.44942 

-.00039 

-.00679 

-.02638 

-.06459 

-.12302 

M5 

-.80036 

-.00441 

.00966 

.08073 

.24023 

.50539 

T5 
x7 

. 0*368 

-.0034C 

-.01372 

-.03040 

-.05176 

-.07367 

■pS 
7 

-.35916 

-.01067 

-.04222 

-.09213 

-.15577 

-.22846 

M^ 

2.60153 

.00944 

.06566 

.20502 

.45221 

.8452? 

x7 

.09909 

.00500 

.01418 

.01910 

.01297 

-.00980 

p8 
8 

.35916 

.01067 

.04222 

.09213 

.15577 

.22846 

Kg 

4.54996 

.19966 

.76128 

1.59939 

2.59827 

3.62984 

n.8 
'8 

.24357 

.02854 

.10847 

.22676 

.36544 

.50081 
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T^"lr  17.  FINAL  ^FWiFR  rMn  STRESSF^ 

UNIT  ITAn  AT  POINTS  74,  75 , 7 . 26 ,77,78 


UNIT  LOAD  AT  POINT 

74 

25 

7 

26 

27 

28 

>i 

.18465 

.15169 

.10918 

.06764 

.03247 

.00860 

Mi 

-7.39595 

-7.07715 

-1.57670 

-1.02849 

-.51705 

-.14269 

Tl 

.31731 

.27277 

.20621 

.13299 

.06604 

.01801 

Fl 

-.18465 

-.15169 

-.10918 

-.06764 

-.03247 

-.00860 

Ml 

-1 .66060 

-1 .76234 

-.82518 

-.46078 

-.19855 

-.04699 

Tl 

-.18765 

-.12997 

-.07370 

-.03288 

-.00988 

-.00113 

*S 

.18465 

.1230? 

.06459 

.02638 

.00679 

.00039 

H| 

.49559 

.35574 

.21100 

.10391 

.03790 

.00713 

T4 
A2 

-.08580 

-.07367 

-.05176 

-.03040 

-.01372 

-.00340 

T?2 
P2 

-0.00000 

.02867 

.04459 

.04126 

.02568 

.00821 

m| 

1.74640 

1.33601 

.87694 

.49119 

.21227 

.05038 

T2 
i2 

-.30794 

-.22577 

-.13821 

-.07103 

-.02803 

-.00599 

F3? 

0.0001  3 

-.02867 

-.04459 

-.04126 

-.02568 

-.00821 

«| 

-1 .74640 

-] .95694 

-1.84227 

-1.38427 

-.76811 

-.22815 

Tl, 

-.30794 

-.35487 

-.34126 

-.25966 

-. 14484 

-.04312 

F5 

.18466 

.22558 

.27554 

.17681 

.10089 

.03057 

m! 

.49^60 

.58197 

.56892 

.43693 

.24462 

.07293 

T5 

.08580 

.07639 

.03504 

-.00520 

-.01391 

-.00640 

»9 

-.18465 

-. 19690 

-.18094 

-.13555 

-.07521 

-.02236 

m2 

1 .66060 

1 .88056 

1.80774 

1.38947 

.78202 

.23455 

Ti 

-.1B765 

-.72709 

-.27766 

-.17726 

-.09978 

-.02981 

*I 

.18465 

.19690 

.18094 

.13555 

.07521 

.02236 

M3 

4 

2.39595 

2.44246 

2.16365 

1.58513 

.86871 

.25632 

'1 

.31731 

.32617 

.29050 

.21176 



.11506 

.03365 
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c 

"MTjiwiJaT  I* 

N  ?F    TAP|  F 

1? 

mm  it  ir:  An 

AT  POINT 

2  4 

75 

7 

26 

27 

28 

P6 

.3153* 

.22846 

.15577 

.09213 

.04222 

.01067 

Vi\ 

-4.  60897 

-3.62984- 

-2.59827 

-1.59937 

-.76128 

-.19970 

m   6 

l5 

.60347 

.50081 

.36544 

.22676 

.10847 

.02855 

F6 

-.31534 

-.22846 

-.15577 

-.09213 

-.04222 

-.01067 

h| 

-1 .5738? 

-.84527 

-.45221 

-.20501 

-.06566 

-.00945 

T6 

-.06829 

-.00981 

.01297 

.01911 

.01418 

.00500 

»« 

-.18465 

-.12302 

-.06459 

-.02638 

-.00679 

-.00039 

Mt 

.7°699 

.50537 

.24022 

.08072 

.00966 

-.00440 

T4 

."8580 

.07367 

.05176 

.03040 

.01372 

.00340 

»z 

B roooo 

.34148 

.22036 

.11851 

.04901 

.01106 

h7 

b 

1 .4880? 

.77161 

.40045 

.17461 

.0^194 

.00605 

^ 

-.72870 

-.49557 

-.25319 

-.09983 

-.02384 

-.00060 

P7 

.50000 

.64852 

-.22036 

-.11851 

-.04901 

-.01106 

h2 

-1  .48802 

-2.76267 

-4.73341 

-2.49487 

-1 .00677 

-.22030 

>; 

-.72870 

-.80568 

-.65204 

-.37087 

-.16284 

-.03932 

»? 

-.18466 

-.22558 

-.22554 

-.17681 

-.10089 

-.03057 

M^ 

.796"9 

.99706 

1 .00983 

.80076 

.46162 

.14103 

>? 

-.08580 

-.07639 

-.03504 

.00520 

.01391 

.00640 

F^ 

-.31  "4 

-.42294 

.44590 

.29532 

.14990 

.04164 

M^ 

1 .^7382 

2.83906 

4.76845 

2.48967 

.99286 

.21390 

rp8 

r7 

-.06829 

-.18138 

-.35779 

-.42989 

-.29878 

-.10177 

F8 

.31 534 

.42294 

.55410 

.70468 

.85010 

.95836 

N§ 

4.60897 

5.46658 

6.13497 

6.46170 

5.90906 

3.90071 

T§ 

.60347 

.64469 

.59874 

.44212 

.23299 

.06379 
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Tr"l  f    }■>, 


tTRFSSPS    at 

in  IT    LOAD    A 


A    <"rrT10N 
T    POINT*    9 


nFTWFN    J 
.10,11.7,1 


UNI' 


UNIT  LOAD 

AT  POINT 

o 

10 

11 

2 

12 

13 

F9 

.07745 

.89568 

.78215 

.64993 

.51628 

.38767 

Hg 

-.67~7'> 

-'.85049 

-5.04264 

-6.09368 

-6.28494 

-5.77115 

T9 

.00071 

-.27636 

-.26390 

-.18938 

-.12620 

-.07123 

?10 

.02745 

.89568 

.78215 

.64993 

.51628 

.38767 

M10 

-.52378 

2.15413 

-.66069 

-2.44785 

-3.38545 

-3.59199 

T10 

.03856 

-.30677 

-.51292 

-.56232 

-.54842 

-.48004 

hi 

.02745 

.10432 

.78215 

.64993 

.51628 

.38767 

MU 

-.37284 

-1 .58828 

3.72629 

1.21662 

-.46018 

-1.38550 

Tll 

-.00O5Q 

.07858 

-.37907 

-.61607 

-.71632 

-.69736 

F12 

.01522 

-.06162 

-.13810 

-.24104 

.36484 

.50000 

M12 

-.13607 

.66538 

1 .73882 

3.45631 

-5.88207 

3.48206 

Tl2 

-.02400 

.06235 

.06356 

-.01837 

.21775 

-.45875 

F13 

.01522 

-.06167 

-.13810 

-.24104 

.36484 

.50000 

H13 

-.04900 

.31551 

.96041 

2.10735 

-3.85435 

6.28304 

T13 

-.03217 

.10517 

.18142 

.22454 

-.20735 

-.03240 

P14 

.01522 

-.06162 

-.13810 

-.24104 

.36484 

.50000 

M14 

.03843 

-.03675 

.17470 

.74235 

-1.79729 

-3.48206 

T14 

-.03263 

.11734 

.23098 

.34896 

-.45410 

-.45875 

P15 

.00057 

-.04055 

-.09471 

-.17161 

-.27048 

-.38767 

M15 

.17119 

-.58716 

-1.08517 

-1 .50048 

-1.67134 

-1.38552 
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c 

NT I  MIOTIC 

H    CF  TABI.F 

13 

UNIT  LOAD 

AT  PCINT 

0 

10 

11 

2 

12 

13 

T15 

-.03661 

.13550 

.27601 

.43481 

.58901 

.69736 

16 

.O0957 

-.04055 

-.09471 

-.17161 

-.27048 

-.38767 

'•'16 

.2268' 

-.8217? 

-1 .63060 

-2.48484 

-3.21707 

-3.59202 

T16 

-.M°?4 

.07399 

.15744 

.26081 

.37558 

.48004 

P17 

.00957 

-.04055 

-.09471 

-.17161 

-.27048 

-.38767 

M17 

.28074 

-1 .05^04 

-2.16362 

-3.45030 

-4.73832 

-5.77118 

T17 

.00292 

-.00773 

-.00822 

.00168 

.02824 

.07123 

P18 

-.01223 

-.04269 

-.07975 

-.10904 

-.11888 

-.11233 

M18 

.06078 

.19511 

.331 16 

.40932 

.39925 

.31097 

T18 

-.00299 

.00009 

.02681 

.09277 

.15780 

.17696 

F19 

-.00565 

-.02107 

-.04339 

-.06942 

-.09436 

-.11233 

Ml  9 

-.0O?85 

-.00385 

.01715 

.07457 

.18400 

.31098 

T19 

-.01124 

0.00000 

-.08081 

-.12354 

-.16056 

-.17696 

P20 

-."l 234 

.04271 

.07911 

.10757 

.11 772 

.11233 

M20 

.22884 

-.79176 

-1  .47569' 

-2.05069 

-7.33232 

-2.30996 

T20 

.03515 

-.11801 

-.21182 

-.27997 

-.29885 

-.27315 

P21 

-.01234 

.04271 

.07911 

.10757 

.1 1772 

.11233 

«21 

.16397 

-.56754 

-1.06093 

-1.48726 

-1.71605 

-1.72263 

T21 

.05230 

-.17736 

-.32257 

-.43444 

-.47561 

-.44921 

P22 

-.01234 

.04271 

.07911 

.10757 

.11772 

.11233 
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CONTINUATION    OF    TARl  F    13 


UNIT  LOAD 

AT  POINT 

Q 

10 

1  1 

2 

12 

13 

H22 

.09785 

-.33900 

-.63810 

-.91251 

-1.08671 

-1.12219 

T22 

.06374 

-.21694 

-.39675 

-.53921 

-.59798 

-.57342 

P23 

-.00011 

.00001 

-.00064 

-.23799 

-.00115 

0.00000 

H23 

.0*622 

-.11684 

-.20142 

-.97848 

-.30890 

-.34016 

T23 

-.01115 

.09825 

.15967 

-.0*556 

.12407 

.02977 

F24 

-.0001 1 

.00001 

-.00064 

-.23799 

-.00115 

0.00000 

H24 

.0?R?6 

-.12491 

-.21773 

-2.14695 

-.32421 

-.34146 

T24 

-.02790 

.08769 

.14137 

-.17202 

.09643 

0.00000 

p25 

-.00011 

.00001 

-.00064 

-.23799 

-.00115 

0.00000 

M25 

.04000 

-.13203 

-.23239 

-3.29908 

-.33706 

-.34016 

T25 

-.02449 

.07648 

.12172 

-.40980 

.06755 

-.02975 

F26 

.00554 

-.02106 

-.04403 

-.07089 

-.09551 

-.11233 

M26 

.08301 

-.29383 

-.56580 

-.82778 

-1.02577 

-1.12221 

T26 

-.03094 

.11400 

.23094 

.36121 

.48495 

.57341 

F?7 

.00554 

-.02106 

-.04403 

-.07089 

-.09551 

-.11233 

"27 

.11277 

-.40665 

-.80121 

-1 .20619 

-1 .53580 

-1.72265 

T27 

-.02239 

.08341 

.1712* 

.27241 

.37311 

.44921 

p28 

. 005^4 

-.02106 

-.04403 

-.07089 

-.09551 

-.1 1233 

M28 

.14165 

-.51637 

-1 .03053 

-1 .57543 

-2.03414 

-2.30998 

T28 

-.01128 

.04311 

.09128 

.15096 

.21724 

.27314 
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•ABIF    14.    CTRFSSFS    AT    A    PFfTirN    BFTWFEN    J  C  T  r>  T r. 
MWIT    LOAD    AT    POINT*    14»3,l*,16.17tlR 


LIMIT  LOAD 

AT  POINT 

14 

3 

15 

16 

17 

18 

p9 

,27048 

.17161 

.09470 

.04055 

.00957 

.44215 

Mg 

-4.73826 

-T. 45032 

-2.16357 

-1.05000 

-.28080 

-3.59172 

T9 

-.02824 

-.00167 

.00821 

.00772 

.00293 

-.02780 

P10 

.27048 

.17161 

.09470 

.04055 

.00957 

.44215 

M10 

-3.21701 

-2.48487 

-1.63060 

-.82168 

-.22689 

-1.12236 

T10 

-.37557 

-.26081 

-.15744 

-.07400 

-.01924 

-.23363 

*n 

.27048 

.17161 

.09470 

.04055 

.00957 

.44215 

Mil 

-1  .67127 

-1.50050 

-1 .08523 

-.58711 

-.17125 

1.35555 

Til 

-.■^BPOO 

-.43482 

-.27602 

-.13551 

-.03662 

-.22345 

P12 

.36484 

.24104 

.13810 

.06162 

.01522 

-.12377 

"12 

1 .70736 

.74243 

.17472 

-.03682 

-.03844 

2.27356 

*12 

-.45407 

-.34897 

-.23097 

-.11733 

-.03264 

-.08672 

p13 

.36484 

.24104 

.13810 

.06162 

.01522 

-.12377 

»13 

3.85441 

2.10743 

.96044 

.31544 

.04900 

1.58573 

T13 

-.20731 

-.22454 

-.18141 

-.10517 

-.03218 

.08178 

p14 

.36484 

.24104 

.13810 

.06162 

.01522 

-.12377 

M14 

5.88213 

3.45640 

1 .73884 

.66531 

.13606 

.88583 

T14 

.21780 

.01838 

-.06356 

-.06235 

-.02410 

.18969 

F15 

-.51628 

.35007 

.78215 

.10432 

.02745 

-.13421 

M15 

-.46002 

6.76504 

-3.72613 

1.58828 

.37289 

-.55746 
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CfNTlNUATICN    CF    TAPI  F    1  4 


UNIT  LOAD 

AT  POINT 

14 

3 

15 

16 

17 

18 

T15 

.71679 

.85833 

-.37913 

.07858 

-.00056 

.17388 

F16 

-.M6?R 

.15007 

.78215 

.1043? 

.07745 

-.13421 

"16 

-3.38*78 

8.60685 

.66085 

7.15420 

.52383 

-1.31516 

Il6 

.^4840 

1.57949 

-.51796 

.24198 

.03859 

.09712 

P17 

-. ■51628 

.35007 

.78215 

.89568 

.02745 

-.13421 

H17 

-6.78478 

10.38317 

5.04281 

7.85048 

.67078 

-2.06285 

T17 

.12620 

2.35861 

-.26393 

-.27636 

.09075 

-.05537 

p18 

-.09436 

-.0694? 

-.04339 

-.02108 

-.00565 

.56592 

Ml8 

.18398 

.07457 

.01714 

-.00381 

-.00385 

-.78853 

T18 

.16056 

.12354 

.08081 

.04070 

.01124 

.00405 

F19 

-.11888 

-.1C903 

-.07975 

-.04269 

-.01224 

.01044 

Mlg 

.39973 

.40929 

.331 19 

.19511 

.06079 

.01036 

T19 

-.15779 

-.09277 

-.02680 

-.00009 

.00298 

-.02532 

P20 

.09551 

.07089 

.04403 

.02106 

.00554 

.30251 

M20 

-7.0340° 

-1.57535 

-1 .0305? 

-.51640 

-.14169 

-2.69794 

T20 

-.21723 

-.15097 

-.09178 

-.0431  1 

-.01 129 

-.02226 

P21 

.09551 

.07089 

.04403 

.02106 

.00554 

.30251 

M21 

-1.53575 

-1.20612 

-.80120 

-.40668 

-.11281 

-1.19181 

T21 

-.37309 

-.27241 

-.17125 

-.08342 

-.02240 

-.19209 

p22 

.09551 

.07089 

.04403 

.02106 

.00554 

.30251 
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c 

MTINUATICN  GF  TAfUF 

14 

UNIT  LCAD 

AT  POINT 

14 

3 

15 

16 

17 

18 

M22 

-1.02572 

-.82770 

-.56579 

-.29386 

-.08306 

.32339 

T22 

-.48493 

-.36121 

-.23094 

-.11400 

-.03095 

-.23001 

p23 

.00115 

.00146 

.00064 

-.00001 

-.00011 

-.13157 

M23 

-.53706 

-.30245 

-.23239 

-.13196 

-.04001 

2.07953 

T23 

-.06755 

-.12629 

-.12173 

-.07649 

-.02450 

-.07063 

P24 

.pons 

.00146 

.00064 

-.00001 

-.00011 

-.13157 

W24 

-.^2421 

-.28309 

-.21774 

-.12485 

-.03827 

1 .42803 

T24 

-.09642 

-.15185 

-.14138 

-.08770 

-.02791 

.08252 

p25 

.00115 

.0C146 

.00064 

-.00001 

-.00011 

-.13157 

M25 

-.30889 

-.26156 

-.20143 

-.11678 

-.03623 

.76565 

T25 

-.12406 

-.17563 

-.15968 

-.09825 

-.03117 

.17829 

F26 

-.11772 

-.10757 

-.07911 

-.04271 

-.01235 

-.12113 

M26 

-1.08660 

-.91246 

-.63810 

-.33908 

-.09789 

-.54281 

T26 

.=■9796 

.53921 

.39675 

.21695 

.06374 

.16926 

P27 

-.11772 

-.10757 

-.07911 

-.04271 

-.01235 

-.12113 

M27 

-1 .71 «Q3 

-1 .48721 

-1 .06093 

-.56762 

-.16406 

-1 .15369 

T27 

. 47*60 

.43444 

.32257 

.17736 

.05231 

.09519 

p28 

-.11772 

-.10757 

-.07911 

-.04271 

-.01235 

-.12113 

M28 

-2.33221 

-2.05064 

-1 .47568 

-.79183 

-.22898 

-1.75578 

T28 

.29886 

.27997 

.21182 

.11801 

.03515 

-.03184 
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T^ULT    t«.    STRFSSF?    at    a    TfTrCN    PFTWFFN    JOINTS 
UNIT    LOAD    AT    POINT*    1 q , 20 , 2 1 . 22 .6 . 23 


UNIT  LOAD 

AT  POINT 

JO 

20 

21 

22 

6 

23 

Pg 

,1465q 

.02236 

.07521 

.13555 

.18094 

.19691 

Hg 

-2.37036 

-.13419 

-.45814 

-.84541 

-1.17680 

-1.36911 

T9 

.03742 

.01659 

.05713 

. 10564 

.14465 

.15976 

P10 

.14659 

.02236 

.07521 

.13555 

.18094 

.19691 

K10 

-1.55124 

-.01106 

-.04407 

-.09930 

-.18098 

-.28527 

*10 

-.13380 

.01025 

.03520 

.06440 

.08537 

.08753 

Fll 

.14659 

.02236 

.07521 

.13555 

.18094 

.19691 

«11 

-.72032 

.11215 

.37033 

.64756 

.81621 

.80075 

Tll 

-.23298 

.01466 

.04945 

.08833 

.11311 

.11004 

P12 

.16186 

-.00821 

-.02568 

-.04126 

-.04459 

-.02867 

M12 

.OR96R 

.18547 

.63533 

1 .17267 

1  .61760 

1  .82141 

T12 

-.20973 

-.02505 

-.08356 

-.14802 

-.19020 

-.18990 

F13 

.16186 

-.00821 

-.02568 

-.04126 

-.04459 

-.02867 

H13 

1 .90227 

.14139 

.49770 

.95218 

1 .38062 

1.67195 

T13 

-.08347 

-.01078 

-.03410 

-.05525 

-.05929 

-.03738 

F14 

.16186 

-.00821 

-.02568 

-.04126 

-.04459 

-.02867 

M14 

2.80039 

.09624 

.35630 

.72444 

1.13312 

1.50977 

T14 

.12185 

-.00041 

.00319 

.01795 

.05046 

.10154 

p15 

-.38872 

-.00860 

-.03247 

-.06764 

-.10918 

-.15169 

"15 

.67724 

-.00083 

.01850 

.08661 

.22267 

.42721 

105 


CONTINUATION    OF    TARLF    15 


UNIT  LOAD 

AT  POINT 

10 

20 

21 

22 

6 

7% 

T15 

.270*2 

•C008R 

-.00040 

-.00899 

-.02795 

-.05620 

Pl6 

-.->R872 

-.00860 

-.03247 

-.06764 

-.10918 

-.15169 

«16 

-1  .50651 

-.04R6? 

-.16170 

-.28824 

-.38153 

-.41117 

T16 

.24^31 

-.00127 

-.00666 

-.01779 

-.03489 

-.05550 

P17 

-.38872 

-.CC860 

-.03247 

-.06764 

-.10918 

-.15169 

H17 

-3.67880 

-.09604 

-.34066 

-.66089 

-.98282 

-1.24642 

T17 

.01691 

-.00759 

-.02859 

-.05923 

-.09446 

-.12787 

p18 

-.01527 

.03057 

.10089 

.17681 

.22553 

.22558 

M18 

.05006 

-.03408 

-.10851 

-.18192 

-.22044 

-.20756 

?18 

.04646 

.00640 

.01391 

.00520 

-.03504 

-.07639 

p19 

.55058 

.00039 

.00679 

.02638 

.06459 

.12302 

Mlg 

-.77262 

.00576 

.01415 

.01159 

-.01416 

-.07482 

T19 

-.05368 

.00340 

.01372 

.03040 

.05176 

.07367 

P20 

.1055? 

.04163 

.85010 

.70468 

.55410 

.42294 

M20 

-1 .67060 

-.84339 

-1 .70873 

-2.99555 

-?. 42744 

-3.41331 

T20 

-.rni62 

.03741 

-.09961 

.02923 

.18124 

.25699 

P21 

.10553 

.04163 

.85010 

.70468 

.55410 

.42294 

Ha 

-1.14208 

-.63785 

2.50460 

.49330 

-.69382 

-1.33407 

^21 

-.13443 

-.02726 

-.06486 

-.08003 

.00130 

.04971 

p22 

.10553 

.04163 

.14990 

.70468 

.55410 

.42294 
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c 

■>NTlwUATir: 

«  OF  TARI.F 

15 

UNIT  LCAD 

AT  POINT 

19 

2U 

21 

22 

6 

23 

M22 

-.60487 

-.42746 

-1.75539 

3.97840 

2.04509 

.75533 

T?2 

-.21070 

-.07377 

-.17879 

.11521 

.06030 

•02444 

P23 

.09026 

-.01106 

-.04901 

-.11851 

-.22036 

.35148 

M23 

.44110 

.16819 

.77510 

1.93256 

3.68400 

-6.02804 

T23 

-.11781 

-.02234 

-.08504 

-.17757 

-.28453 

.39960 

P24 

."9026 

-.01106 

-.04901 

-.11851 

-.22036 

.35148 

M?4 

.R9543 

.11488 

.53753 

1 .35543 

2.60655 

-4.30417 

T24 

-.05945 

-.00998 

-.02773 

-.03401 

-.00988 

-.05151 

P25 

.09026 

-.01106 

-.04901 

-.11851 

-.22036 

.35148 

H25 

1 .34294 

.06069 

.29587 

.76798 

1 .50926 

-2.54755 

T25 

.03828 

-.00232 

.00866 

.05870 

.16982 

-.35066 

^26 

-.35916 

-.01067 

-.04222 

-.09213 

-.15577 

-.22846 

"26 

.80923 

-.04372 

-.14433 

-.25240 

-.31990 

-.28537 

*26 

.24891 

.00350 

.01075 

.01703 

.01875 

.01464 

P27 

-.35916 

-.01067 

-.04222 

-.09213 

-.15577 

-.22846 

M27 

-.98922 

-.09656 

-.35321 

-.70790 

-1.08957 

-1.41379 

T?7 

.24105 

-.00262 

-.01098 

-.02490 

-.04279 

-.05954 

p28 

-.'K"16 

-.01067 

-.04222 

-.09213 

-.15577 

-.22846 

M28 

-2.78015 

-.14865 

-.55941 

-1 .15801 

-1 .85095 

-2. 51145 

T28 

.07647 

-.01333 

-.05082 

-.10636 

-.17118 

-.23180 
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"API  F  16.  STRESSES  AT 
UNIT  LCAO  A 


A  crcTICN  RFTWFFN  JOINTS 
T  P~INT<;  74, 7^, 7 .76.77.78 


UNIT  LHAD 

AT  PCINT 

7* 

75 

7 

76 

27 

28 

r9 

.1846* 

.15169 

.10918 

.06764 

.03247 

.00860 

Mg 

-1  .38996 

-1 .24639 

-.98289 

-.66087 

-.34068 

-.09600 

T9 

.15201 

.17788 

.09446 

.05923 

.02859 

.00759 

F10 

.18465 

. 15169 

.10918 

.06764 

.03247 

.00860 

"10 

-.37338 

-.41113 

-.38159 

-.28821 

-.16171 

-.04858 

TlO 

.07503 

.05551 

.03488 

.01779 

.00666 

.00128 

*U 

.18465 

.15169 

.10918 

.06764 

.03247 

.00860 

Mll 

.64603 

.42724 

.22261 

.08663 

.01848 

-.00079 

Tll 

.08693 

.05621 

.02794 

.00899 

.00040 

-.00088 

P12 

o.onoon 

.02867 

.04459 

.04126 

.02568 

.00821 

M12 

1  .76650 

1 .50968 

1 .13306 

.77444 

.35639 

.09626 

l'l2 

-.15456 

-.10157 

-.05045 

-.01795 

-.00320 

.00041 

F13 

0.00000 

.02367 

.04459 

.04126 

.02568 

.00821 

M13 

1 .77334 

1.67186 

1.38055 

.95218 

.49780 

.14141 

T13 

-0.00000 

.03738 

.05929 

.05525 

.03409 

.01079 

P14 

P. 00000 

.02867 

.04459 

.04126 

.02568 

.00821 

"14 

1.76659 

1.82131 

1 .61754 

1.  17268 

.63542 

.18549 

T14 

.15455 

.18990 

.19019 

.14802 

.08357 

.02506 

p15 

-.18465 

-.19690 

-.18094 

-.13555 

-.07521 

-.02236 

M15 

.64610 

.80069 

.81626 

.64753 

.37044 

.11219 

108 


c 

r-NTj  nuatio 

N  CF  TARLF 

16 

UNIT  LOAD 

AT  POINT 

?4 

25 

7 

26 

27 

78 

T15 

-.08694 

-.1 1001 

-.11317 

-.08832 

-.04946 

-.01467 

P16 

«.]8465 

-. 19690 

-. 18094 

-.13555 

-.07521 

-.07236 

M16 

-.373^1 

-.28526 

-.18094 

-.09933 

-.04396 

-.01102 

T16 

-.07503 

-.08752 

-.08538 

-.06439 

-.03521 

-.01025 

p17 

-.18465 

-.19690 

-.18094 

-.13555 

-.07521 

-.02236 

M1? 

-1.38988 

-1.36905 

-1.17675 

-.84544 

-.45803 

-.13415 

T1T 

-.15201 

-.15975 

-.14466 

-.10564 

-.05713 

-.01659 

P18 

.18465 

.12302 

.06459 

.02638 

.00679 

.00039 

M18 

-.15069 

-.07483 

-.01417 

.01158 

.01414 

.00577 

he 

-.08580 

-.07367 

-.05176 

-.03040 

-.01372 

-.00340 

P19 

.1.8466 

.22558 

.22554 

.17681 

.10089 

.03057 

M19 

-.15071 

-.20756 

-.22047 

-.18191 

-.10850 

-.03407 

T19 

.08580 

.07639 

.03504 

-.00520 

-.01391 

-.00640 

P20 

.31534 

.27846 

.15577 

.09213 

.04222 

.01067 

M20 

-3.08675 

-2.53144 

-1 .85098 

-1.15807 

-.55934 

-.14874 

T20 

.26747 

.23180 

.17118 

.10637 

.05081 

.01334 

P21 

.31534 

.22846 

.15577 

.09213 

.04222 

.01067 

«2i 

-1.54103 

-1.41378 

-1.08959 

-.70796 

-.35314 

-.09664 

T21 

.06542 

.05955 

.04279 

.02490 

.01097 

.00262 

P22 

.31534 

.22846 

.15577 

.09213 

.04222 

.01067 
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CONTINUATION    OF    TARLE    16 


UNIT  LOAD 

AT  POINT 

?4 

25 

7 

26 

27 

?8 

M22 

.01641 

-.28536 

-.31992 

-.25246 

-.14425 

-.04381 

T22 

-.00115 

-.01464 

-.01875 

-.01704 

-.01075 

-.00351 

F23 

.^oooo 

.34148 

.22036 

.11851 

.04901 

.01106 

M23 

4.P150B 

?.49B?4 

1 .509?? 

.76790 

.?95R5 

.06070 

T23 

-.4RRA3 

-.'5280 

-. 16981 

-.05868 

-.00865 

.00?31 

P24 

.^nnoo 

.'4148 

.22036 

.11851 

.0490] 

.01106 

M24 

6.51 ]5fl 

4.20585 

2.60651 

1 .35534 

.53751 

.11488 

T24 

-.0288? 

-.06010 

.00989 

.03402 

.0?773 

.00998 

P25 

.50000 

.34148 

.22036 

.11851 

.04901 

.01106 

M25 

-4.01508 

5.88145 

3.68396 

1.93247 

.77508 

.16820 

T25 

-.48843 

.38032 

.28453 

.17757 

.08504 

.02234 

F26 

-.31534 

-.42294 

.44590 

.70468 

.14990 

.04164 

M26 

•0165C 

.75541 

6.98353 

-3.97832 

1.75539 

.42759 

T26 

.00114 

-.02444 

.15531 

.11518 

-.17879 

-.07376 

P27 

-.31534 

-.42294 

.44590 

.70468 

.14990 

.04164 

M27 

-1  ,c&n<)4 

-1 .33399 

9. 1 4546 

-.493?2 

?.K0456 

.63803 

T27 

-.06541 

-.04970 

.8595? 

-.08005 

.00721 

-.02724 

P28 

-.31534 

-,U->2<iti 

.44590 

.70468 

.85010 

.04164 

M28 

-3.08666 

-3.41323 

11 .23778 

2.99564 

1 .70873 

.84361 

128 

-.26746 

-.25697 

1.74947 

.02921 

-.09961 

.03745 
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Since  a  problem  in  the  field  of  gridded  frameworks  can  be 
reduced  to  the  solution  of  a  set  of  simultaneous  linear  equations, 
the  difficulty  in  solving  these  problems  involves  the  solution  of 
a  set  of  simultaneous  equations.   If  any  of  the  members  of  the 
framework  are  curved  in  plan,  many  trigonometric  terms  will  appear 
in  the  coefficients  of  these  equations.   This  makes  the  analyses 
of  the  problems  in  this  field  very  complicated.   Efforts  to  solve 
this  kind  of  problems  have  been  made  by  many  authors.   Most  of 
them  have  tried  to  use  the  moment  distribution  method  to  avoid 
solving  a  set  of  simultaneous  equations.   However,  the  allowed 
simplifications  will  be  much  reduced  if  the  structure  is  of  the 
gridded  type.   Since  the  moment  distribution  method  itself  will 
become  tedious  if  it  is  used  to  solve  the  problems  of  the  gridded 
type. 

In  this  report,  advantage  will  be  taken  of  electronic  computer 

procedures  to  solve  the  simultaneous  equations  and  to  do  the 
analysis  work.   It  allows  the  greatest  possible  simplification  in 
the  computation  required.   The  member  can  have  any  curvature  in 
the  plan,  including  segmental,  but  particular  attention  will  be 
paid  to  the  common,  circular-arc  shape. 

Formulas  are  presented  for  frames  of  constant  cross  section 
that  give  the  relationship  between  the  end  moments,  end  torques 
and  shears  and  the  corresponding  rotations  and  displacements  of 
the  end  sections.   These  relations  have  been  derived  from  the 
expression  for  the  total  strain  energy  due  to  bending  and  torsion 
by  means  of  Castigliano's  theorem.   From  these  relations  the 
stiffness  matrix  of  each  single  member  can  be  written  and  the 


stiffness  matrix  for  the  whole  structure  follows.   Through  compu- 
tation with  the  electronic  computer  the  end  forces  can  be  obtained. 
These  end  forces  are  combined  with  the  fixed-end  forces  to  ob- 
tained final  end  forces. 


